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ABSTRACT 

A research program into three aspects of space satellite dynamics 
has been carried out. 

First, a four-dimensional space-time formulation of Newtonian 
mechanics is developed. This theory allows a new physical interpretation of 
the conservation theorems of mechanics first derived rigorously by Noether. 
The formulation has turned out to be similar to that in a forgotten 1923 paper 
by E. Cartan. However, the work presented here offers much greater 
physical insight into the underlying mathematical structure of Newtonian 
mechanics than that of Cartan. 

Second, a new concept for estimating the three angles which specify 
the orientation in space of a rigid body is presented. Two separate methods 
for implementing this concept are discussed, one based in direction cosines, 
the other on quaternions. Two examples are discussed: constant orientation 
in space, and constant rate of change of the three angles with time. The 
behavior of each method in the absence of noise is discussed. Fokker-Planck 
equations are derived for the aposteriori probability density functions for 
estimation error in the presence of noise. Steady state irror statistics for 
the case of constant rate of change of the three angles with time for a 
quaternion implementation are explicitly derived. 

Third, two synchronous equatorial orbit communication satellite 
designs which use sunlight pressure to control their attitude are analyzed. 
Each design is equipped with large reflecting surfaces, tailed solar sails, 
which can be canted in different directions to generate torques to correct 
pointing errors. The total equations of motion for each design are derived, 
and then linearized about a nominal trajectory; a theoretical analysis of the 
linearized equations is carried out. A specific control Lew is discussed. 
Disturbance torques are shown to be negligible compared to attitude control 
torques. Computer simulation of the total equations of notion verified the 
theoretical analysis of the linearized equations. 


New directions in space satellite dynamics research are sketched. 
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Chapter One 
Introduction 

1.1. Introductory Remarks 

In 1958 the first communication satellite (SCORE) was placed in earth 
orbit. Since that time, dozens of communication satellites have been launched 
into space and are currently being used for commercial and governmental 
communication between points scattered across the face of the earth. 

One fundamental problem in communication satellite design is control 
of satellite orientation. This report describes the results of a three fold 
research program into space satellite dynamics. The first and second phases 
deal with aspects of deterministic and stochastic rigid body mechanics, 
respectively. The third phase builds on the first two, and involves the analysis 
of two specific satellite designs which use sunlight pressure for attitude control. 

1.2. Scope of the Research Program 

1.2, 1*. Theoretical Rigid Body Mechanics 

The first phase of the research program was devoted to applications of 
differential geometry, tensor algebra and calculus, and Lie groups and Lie 
algebras to problems in space satellite dynamics. The main result was the 
development of a four-dimensional space-time formulation of Newtonian 
mechanics. This theory is a complete self-consistent set of ten equations; 
four of the equations deal with linear momentum, while the other six are con- 
cerned with angular momentum. Many mechanics treatises often ignore three 
of the angular momentum equations; this is not surprising, since information 
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contained in these three equations is also contained in the linear momentum 
equations, making these equations to a certain extent redundant. However, 
by including these three additional equations, a new physical interpretation 
of the ten conservation theorems of mechanics due to Noether arises in a very 
natural simple manner. Noether's approach to the conservation theorems 
was based on a Lagrangian formulation of mechanics plus some subtle varia- 
tional arguments; the same theorems follow quite naturally and simply by 
formulating Newtonian mechanics in space-time. No attempt is made to 
achieve mathematical rigor in developing this theory of mechanics; rather, 
the emphasis is on physical insight and intuitive arguments, coupled with the 
requirement the theory be consistent with experimental observations. 

1.2.2. Estimation Theory 

The second phase of the research program touched on a problem in 
estimation theory: given noisy measurements of the orientation in space of a 

rigid body, estimate the three angles which define this orientation. 

Is attitude estimation really a problem in space satellite dynamics? 

Yes, in fact in many present day satellites, attitude estimation is the chief 
limitation in controlling the spatial orientation of these designs; the reader 
should consult the bibliography for specific examples. One approach to attitude 
estimation is to linearize the equations of motion about a nominal trajectory, 
and then apply linear filtering theory to the linearized equations of motion. 

The actual implementation of this approach can be very complex. In addition, 
this approach ignores the structure implicit in the equations of motion which 
describe how orientations in space change with time. 
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The method presented here is analogous to a phase-locked loop, in that 
it estimates the three angles defining the spatial orientation of a rigid body 
using an estimation procedure resembling that of a phase-locked loop. This 
method takes advantage of the structure implicit in the equations of motion, 
and for the two examples considered here (constant orientation in space, and 
constant rate of change of the angles with time) is quite simple to implement 
compared to an approach based on linearization of the equations about a nominal 
trajectory. 

The chief question left unanswered is what is the optimal or best method 
for processing noisy sensor measurements in order to estimate spatial orien- 
tation. The method presented here works well for the two examples considered; 
however, it is still not clear how well the optimal estimation procedure would 
perform compared to the method discussed in this report, and what if anything 
is lost using the method de scribed here or a method based on Kalman filtering, 
1.2.3. Solar Pressure Attitude Control 

The attitude control dynamics of two specific satellites which use sun- 
light pressure to generate attitude control torques are analyzed in the last 
phase of the research program. 

Why use sunlight pressure for attitude control? At synchronous altitude 

the largest disturbance torque on many present day communication satellites is 

N, 

due to sunlight pressure (NASA(63) ); since it is such a nuisance, perhaps it can 
be used to aid rather than hinder attitude control. By canting large reflecting 
surfaces, called solar sails here, in specified directions, it appears to be 
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possible to control the orientation of two specific designs with as great an 

) 

accuracy as any other presently existing method. 

All of these statements are based on paper-and-pencil analysis plus 
extensive computer simulation of each design. Clearly, this is no substitute 
for an actual test in space of these ideas. Some attempt was made to include 
technological constraints into the designs de scribed here , and all computer 
simulations were carried out with hopefully realistic numbers for all param- 
eters; however, many engineering problems were ignored, on the grounds 
that they did not critically affect the attitude control dynamics of each design. 

Thus, the research program has merely indicated that solar pressure attitude 
control of synchronous orbit communication satellites might be feasible, but 
the question still deserves more discussion. 

1.3. Background 

The reader is referred to the bibliography for material that was found 
to be especially helpful in one or more aspects of the present research program. 

In particular, Goldstein( 1 0) provides sufficient background in classical 
rigid body mechanics, while Abraham(l) provides a much more modern treat- 
ment of mechanics . Nelson(20), Flanders(7), Greub(ll), and Warner(27) 
provide sufficient background material in tensor algebra, as it is used in the 
first part of this research program. 

Van Trees(47), Viterbi(48) and Jazwinski(39) are excellent references 
in estimation theory and phase -locked-loop techniques. Ito(36) and McKean(4 1 , 44) 
provide a more modern and mathematically rigorous approach to the character- 
ization of noisy measurements of spatial orientation. 



To gain some idea of present day attitude control techniques that are 
used, the reader is referred to Likins(58), Fleischer(55) and Much et al(62). 
The NASA publication on radiation pressure torques, including both sunlight 
pressure torques as well as thermal torques, is a good introduction to the 
practical aspects of solar pressure attitude control disturbances(NASA(63) ) . 
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Chapter Two 

Four -Dimensional Space -Time Newtonian Mechanics 
2 . 1 Introduction 

In this chapter, a novel, simple and elegant space -time formulation 
of Newtonian mechanics is developed. The main mathematical tool used to 
develop this theory is the exterior product, a generalization of the cross 
product between two vectors from 3 to 4 dimensions. Using this tool, a 
complete and self-consistent set of equations of motion for a single particle 
are developed; four equations deal with time derivatives of linear momentum, 
while the other six are time derivatives of angular momentum. The first 
four, Newton's laws, are the fundamental equations of motion; the six angular 
momentum equations follow as a natural consequence of the four linear 
momentum equations. 

The formulation presented here offers a new simple physical interpre- 
tation of the ten conservation theorems of mechanics, theorems dealing with 
ten quantities that do not change with time when no forces are present. These 
theorems were shown by Noether(21) to be the only possible conservation 
theorems; however, her arguments, based^on a Lagrangian formulation of 
mechanics, involve some subtle variational techniques viz. a. viz. the ten 
independent parameters governing transformations from one space-time 
coordinate frame to another. On the other hand, these ten theorems follow 
quite directly and simply from the ten equations to be presented here. 

The theory is complete and self-consistent, in that no more equations 
of motion can be found than those presented here, provided the only 
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mathematical operations allowed are the exterior product and differentiation 

I 

with respect to time. Many treatises on mechanics hav overlooked three 

of the angular momentum equations to be presented here, and concentrated 

- •- 

on the four fundamental linear momentum equations, often leaving the im- 
pression these three angular momentum equations do not even exist. After 
this work had been completed, the author searched the literature and was 
able to find only one forgotten paper by Cartan(3) which even intimated there 
might be ten equations of motion for Newtonian mechanics (see Appendix). 
However, Cartan's work is very different in spirit from that found here. In 
this chapter, a series of intuitive and physical arguments are presented to 
develop a natural and correct approach to Newtonian mechanics, at the 
expense of mathematical rigor. Cartan's treatment contains more mathe- 
matical rigor, but is missing much of the physical flavor found in the argu- 
ments here; moreover, he is extremely terse on why there should be ten 
equations of motion for Newtonian mechanics. Since the three forgotten 
angular momentum equations contain no new information than that found in 
the four linear momentum equations, it is perhaps not surprising they have 
been forgotten. 

It is straightforward to extend these equations to include effects due 
to special relativity in a natural manner, unlike other relativistic theories 


of mechanics known to the author; this issue will be dealt with elsewhere. 
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2.2. Kinematics of a Single Particle 

i 

This section is concerned with describing the position of a single 
particle, in space and time with respect to a reference coordinate. Space 
and time together make up the arena in which the dynamics of the single 
particle can take place. The particle is assumed to have positive mass M, 

i 

but occupies such an infinitesimal volume of space at any instant of time that 
1 it can be considered a point in space-time with mass, or a point mass. 

Four independent numbers are needed to describe the position in space and 
time of the particle: three to specify its position, in space, and one to 
describe the instant of time the particle is at that position. 

[ 

An interesting question now arises: given the space-time coordinates 


of a particle in a reference coordinate frame, called frame A (arbitrarily), 
how do these coordinates relate to the particle's coordinates in a different 
reference frame, labeled frame B (again arbitrarily)? It is assumed both 
frames have a standard orthonormal rectangular right-handed set of basis 
vectors for the three space coordinates; the unit vector associated with time 
is assumed orthogonal to the three space basis vectors. The relationship 
between the particle's coordinates in frames A and B is (Goldstein( 10) ): 
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where (x ,y ,z A ,t ) and (x , y , z t ) are the coordinates of the point 
A A A A B B B B c 

mass in space-time in frames A and B respectively. 

The physical interpretation of the various parameters in this coordinate 
transformation is now discussed: 

*> ,X o- WV-- when ( *A = °- y A = °' 

Z A ‘a 0) ’ then (X B = V y B = V 

z n = z » = t )' where (x ,y ,z , t ) 

r> o B o oooo 

are all assumed constant. Thus, the 

space-time origin of frame B is displaced 

from the space -time origin of frame A by 

(x , y , z , t ) 
o o o o 

ii) (v , v , v ) - - when (x A = 0, y A = 0, z A = 0, 
x y z A A A 

t. = t ) , then (x = x + v t , y = y + 

A A B o x A 7 B 7 o 

V t , z = z + V t , t = t + t ). This 

yBB o zzB A o 

means the spatial part of the space -time 

origin of frame B is displaced from the 

spatial component of the space -time of 

frame A by the sum of (x ,y , z ) and 

o o o 

(v t ,v t. ,v t ). The first term is 
x A y A z B 

part of the displacement discussed in i); 
the second term is called a translation. 


The time components t and t are re- 
lated by a displacement in time, t . 

o 
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Unlike (x , y , z , t ), (v , v , v ) are assumed 
o 7 o o o x y z 

to be functions of time . 


iii) (d 11 ,d 12 ,d 13 ,d 21 ,d 2r d 23 ,d 32 ,d 33 )--li.e S e nine 
parameters have six constraints associated with 
them: 
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constraints 


where d^represents the projection of a unit vector directed along frame A's 
x-axis onto frame B's x-axis, d the projection of a unit vector directed along 

X 

frame A's y-axis onto frame B's x-axis, and so on. Since there are a total of 
nine parameters with six constraints, there are only three independent 
parameters associated the d 's,i,j = 1,2,3. To emphasize that there are 
only three independent parameters, the matrix of d. ,'s i j = 1,2,3 is decom- 
posed into the product of three matrices, each parameterized by one independent 
variable, each representing a rotation in space about some axis: 
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* xp(a lil ) ' Xp(a 2 k> ) Lj+L^L^L, 


0 0 0 

0 0 -1 

0 1 0 


L = 


0 0 1 

0 0 0 

-10 0 


L = 
— z 


0 -1 0 

10 0 
0 0 0 


where Exp(a L ) k = 1,2,3 is defined as a matrix exponential. 

iC K 

Thus, the ten independent variable s , a,,a_,a_,v ,v ,v ,x , y ,z ,t 

12 3xyzo ooo 

completely describe the relationship between space -time coordinates in 
frame B with respect to those of frame A. 

2.3. Exterior Algebra 

In order to discuss Newtonian mechanics in four dimensions, the 
notion of a cross product between two vectors must be extended from three 
to four dimensions. This extension is called an exterior product. The treat- 
ment here will concentrate on the algebraic properties of the exterior 
product; for much more detailed discussions of exterior algebra, the reader 
should consult Flanders( 7 ), Warner(27), Nelson(20) or Greub(ll). Exterior 
products will be defined in very general terms with respect to a finite dimen- 
sional vector space; the special case of a four dimensional vector space, 
space-time, will be used as a concrete non-trivial illustration of these general 
propertie s . 

Let V be a real n-dimensional vector space, with an orthonormal set 
of basis vectors {dx , . . . , dx } with re spe ct to some coordinate frame . 



25, 


Associated with this vector space is a space of zero-dimensional vectors, 

I 

the space of all real scalar functions of the n coordinate j (usually these func- 
tions are assumed to be continuous and infinitely differentiable, but these 
properties will not be needed here). This zero- dimensional vector space is 
denoted AN) , and has a basis vector {l}. Next, there exists an n-dimensional 
vector space, with the same basis as V, denoted A'M . Third, there exists 
an ^ n(n - 1) dimensional vector space associated with V, denoted m 
iisting of all vectors of the foi 


consis 


>rm 


2.4) -vT= ; 


n * 

z. z 


r -« 


r -*j Ux.Ackp ^ & R 


where "a” the exterior product. The exterior product must obey the 
following four constraints: 

t) (C.TT, f A c, KjAu) +c,j. (fj Au) 

li) < A d t (>rAU,) + d a (TAMj 

Hi) v A V = 0 

VAV4 M.AV >0 

^ y\ * 

>T* Z 4*, = <5 V. dx; , V Z ^ 

;»» •-» 

n m ( * j 

,5 u*dx; , = Z u ] ^ Z u j 


u ; ,n V'*i e ^ 


Thus, a basis for A(V) is {dx.Acbc. ; i,j = 1 ,...n} subject to the constraints of 
exterior multiplication; there are ^ n(n-l) vectors in this basis forA.^ 
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By induction, the above rules can be extended to vector spaces /M 
where 2 ^ p ^n. Each vector v in m can be written 


as 


A * 

2.5) V* 5 

'l' » 1-t 0 X , 

Vectors in A*W) are subject to the following conditions 


1) (av + bu) Ar A a- ^A^A— A** p ) + b (<A*“ 4 A~- Ar p ) 

v t A U^+b-u) Av 3 A— Av f -- ^K^,V- A/p) + bKAXTAiTjA - AVpO 
A “ A *p-' A Cav+bu)=. aCv, A -av^Av) + b C^A-Avy* Avi) 

W V', A-Av-,’ 0 IF FOR 50ME MR OP uoDitt&S f.--/. Poft 

ill) v“ t A-> AVp c^AW^S Sl^M (F FOR -SOME PAift Op AD\TAC£dl 

AftE 


where a, V", V*j } — ,Vp £ /'(V) Q t b e R 1 

As an example, consider the four dimensional space -time vector space 
discussed earlier, with an orthonormal rectangular set of basis vectors, 

{dx, dy , dz , dt}. Associated with this vector space are a total of five vector 
spaces, defined according to the exterior multiplication rules of exterior 
algebra: 


Vector Space 

Basis Vectors 

A°(v) 

i 

A 'ey) 

»dt 

A J CV) 

okj/tfj, dyicfc ,dx Ady, dt 

A*CV) 

o/^AchjAcit 


elx A^ AdjAdt 
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The physic- -terpretation of these basis vectors is of interest: 

i) ( l} a scalar real number, with neither 

magnitude nor direction in space -time 

ii) A*tV)- - {dx, dy , dz , dt} are basis vectors with both 
magnitude and direction, e.g. , dx has a "one- 
dimensional magnitude" of+1 in the +x-direction 

iii) A*(V) - {dy A dz , dz A dx, dxAdy , dtA dx, dt Ady, dtA dz } 
have both magnitude and direction, e.g. , dy/dz has 

a "two-dimensional magnitude" of+1 in the "y-wedge-z" 
direction, while dt dx has a "two-dimensional magnitude" 
of +1 in the "t-wedge-x" direction, where "wedge" denotes 
exterior product and not the ordinary three dimensional 
vector cross product 

iv) aW- - { dyA dz A dt , dz A dxA dt , dxA dy A dt , dxA dy A dz } have 
both magnitude and direction, e.g. , dyAdzAdt has a 
^three-dimensional magnitude" of + 1 in the "y-wedge z- 
wedge-t" direction, while dxAdyAdz has a "three-dimen- 
sional magnitude" of+1 in the "x -wedge -y-wedge-z" 
direction 

v) a"cv) -- {dxA dyAdzAdt } has a "four dimensional magnitude" 
of+1 in the "x-wedge-y-wedge-z-wedge-t" direction. 

2.4. Newtonian Space -Time Mechanics of a Single Particle 

This section derives the ten space -time equations of motion of a 
single particle the space -time coordinates of the particle are 
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2 . 6 ) 


= space -time coordinates of a single particle 


The velocity of the particle is defined to be the derivative with respect to 
time of the space -time coordinates: 
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= velocity of a single particle 


The linear momentum of the particle is defined to be its mass in times its 
velocity: 
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= linear momentum of a single 
particle 


Newton postulated that the time rate of change of the linear momentum of a 
single particle equals the forces acting on the particle: 
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Where (F ,F ,F ,0) are forces acting on the x-,y-,z- and t- components of 
x y z 

» 

the linear momentum, respectively, and these equation assume the linear 
momentum and forces are computed in a special c .rdinate frame, an 
inertial frame, which is defined (albeit circularly, as a coordinate frame in 
which the linear momentum equations of motion can be written in the form 
above . 

The moment about the space -time origin of the linear momentum can 
be found using the exterior product: 

+ A 4 M<itD = 

2.10) ^Adx i ~'$0pv) f 

'WvXo^ olt 4 Ctgpjj-*k^<0 4 Adj. 

The first three terms, along dyAdz,dzAdx and dxAdy, are called here space- 
like angular momentum components, because they arise from purely spatial 
exterior products. The final three terms, along dt Adx, dtA dy , dtAdz are 
called here time-like angular momentum components, because they arise 
from a mixture of space-time exterior products. All six terms together 
make up what is called here the total angular momentum of the particle. 

Using the four equations of motion for linear momentum, the time 
rates of change of the six angular momentum components are: 


2.11) 
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t 0 F x + (*<*)- C*«a x ) ~(^) x 0 2 -L 0 Fn 
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2.15) t,F } 

2.16) &(V>-"P = t.Fjt t„Fj 

Next, all possible exterior products of (x dx + y dy + z dz + t dt) 

o o o o 

and (p dx + p dy + p dz + m dt) taken three at a time are computed, and 
x y z 

all eight such products are found to be zero. Thus, there are no further 
quantities than the four components of linear momentum and six components 
of angular momentum in mechanics, given all that is computed are exterior 
products of space-time vectors describing position and linear momentum, 
and their time derivatives. 

What if no forces act on the particle? IfF = F = F = 0, then all 

x y z 

components of the linear and angular momentum are constant. Note that 

while it is obviou^ from these ten equations that quantities such as p or 

(y p - z p ) are now constant, it is also obvious that (t p - mx ) is 
o z o y oxo 

constant; the first two constants appear in many texts on mechanics, but the 
last one is mentioned rarely. 

Noether(21) derived these ten conserved quantities from the Lagrangian 

L. of the particle, assuming no forces act on the particle, where 

1 2 2 2 

L = — m(u 4- u + u ) 

2 x y z 

through an argument based on techniques used in calculus of variations. 

Recall that four number specify the position of a particle in space -time; if 
this position is perturbed slightly, the new position will be a function of the 
old position as well as the ten independent parameters defining the transfor- 
mation from the old to the new position. This transformation is 
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Associated with each independent parameter, according to Noether, there 
is a conserved quantity 

Parameter . Conserved Quantity 
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D*Y s»b 
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t>oPx - XoP^ 

A<t^ 

Mo?# 

AV* 

t#Px ~*\Xo 

AVj 





where E is the kinetic energy of the particle, 

EL = L - £ Pi + 1 ~ 

Note that since p ,p ,p and E are constant, m must be constant. Note 
x y z 

further that E equals the Lagrangian L; Noether demanded that the ten inde- 
pendent parameters not perturb the Lagrangian, but due to the fact that the 
kinetic energy and Lagrangian are identical here, perhaps missed a more 
fundamental observation: m is constant, a fact that arises quite naturally 
from a space-time formulation of mechanics. 
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2.5. Space -Time Equations of Motion of a Single Rigid Body 

This chapter concludes with a derivation of the ten space-time equa- 
tions of motion for a single rigid body. A single rigid body may be viewed 
as finite collection of small bodies, each occupying a small but finite volume 
of space at any instant of time, and fixed rigidly with respect to all the other 
small bodies comprising the rigid body. "Small" in this context means the 
volume occupied by each body comprises many thousands of atoms of what- 
ever substance the rigid body is made of, while at the same time this volume 
is much less than total volume of the rigid body. The sketch below shows 
a typical rigid body decomposed into N small bodies. 

Figure 2.1 

A single rigid body 
broken down into 
N rigid bodies 





th 


The space -time coordinates, in an inertial frame, of the k— mass are 
(x^,y^,z^, t^). The equations of motion for this mass, assuming it is a 
point mass, are: 
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2 . 20 ) 

2 . 21 ) 

2 . 22 ) 
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2 ! 27) 
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F^ , F ^ , F ^ = x- , y-, and z- components , respectively, of 

xk yk zk 

constraint or reaction force of mass j acting 
mass k 

if mass k is adjacent to mass j 
if mass k is not adjacent to mass j 


kj 


f ' 

10 i 


If all equations for the x- component of the linear momentum are added. 


JJ Jf 


defining 


i ( <r p > - ( 2 f ) + ( 2 2 f j s:J 
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2.30) F j <5 F 




X* 


and realizing the reaction force of mass k on mass j is equal and opposite to 
the reaction force of mass j on mass k, 
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2.32) 2 ^ " 0 
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This equation can be written simply as 


d „E 

2.33) — P = F 

dt x x 

By adding up all the other equations for the linear momentum, component by 
components, it is easy to see 
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To proceed further, it is useful to define the center-of-mass of the body. 
This is a set of coordinates, (x ,y , z ,t ), where 
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Note that time is assumed identical through the body, at all space-time 
coordinates. It is also useful to define coordinates with respect to the center 
of mass: 

2 . 41 ) 

2.42) ‘jilt* ' (Jf 
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Then in terms of center-of-mass coordinates, the linear momentum equa- 


tions of motion become: 
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Adding up the angular momentum equations of motion, component by 
component, it is straightforward to show 

2 • 4 9) jfc ( ‘t'oM. - A X CM^ - ^0** 

2.50) fy” A ^CM 

2.51) j*. (tcM F J 

j j JX ^ 

2.52) *{ ^ j = 2 1 > P J|C 3 

2.53) ^5 £ 1 = £ E)f F X« - X « F yJ 

L £3 1 w K 3 * 

2.5 4 ) £ j £ I *« p.,- Pkk 3 j = £ [ *« F i, - > F * ( . 3 

1 £31 • |c-»l 9 

These ten equations, (2.45)*(2 . 54) are the ten space-time equations of motion 
for a single rigid body. Note that the time -like component of the angular 
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momentum is a function of center-of-mass variables alone, as is the linear 

1 

momentum, while the space-like component of the anguxar momentum depends 
both on center-of-mass variables and coordinates with respect to the center 
of mass . 

The equations of motion decouple into two sets of equations, one 
describing motion of the center of mass (equations (2 .45)-2 . 5 1 ) ) , the other 
motion about the center of mass (equations (2.52), (2.53), (2.54)). Since 
the motion of the center of mass in space can be described by six first order 
differential equations, three of the equations (in (2.45) - (2.51) ) are redundant; 
typically, equations (2.49), (2.50), (2.51) are ignored, since they contain no 
more information about how external forces act on the center of mass than 
equations (2.45a), (2.46a), (2.47a). In fact, the motion of the center of mass 
is completely described by 
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APPENDIX 


After the work described in Chapter II was completed, a thorough 

literature search was undertaken to see whether a similar four -dimensional 

space-time formulation of Newtonian mechanics had been presented elsewhere. 

To the best of the author's knowledge, there is only one earlier work that is 

similar to the theory presented in Chapter II: 

M Sur les Variety s k Connexion Affine et la Theorie de la 
Relativite Gen£ralisee ," E . Cartan, Annales 
Scientifique s de l'Ecole Normale Superieure, Series 
3, Volume 40, 1923 

Since this paper has apparently been forgotten, the three pages of the article 
relevant to the work described in Chapter II are included here so that the 
reader may compare the two developments for himself. 
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dies sont idcutiquement veriliees s’il u’y a pas de pression; dans le 
cas general dies donnenl 

Pz- — 1>SZ— ••• 

Pxz—Pz^—O, 

P,* — P*, = 0. 

II. On peut re presenter les resuifats precedents an moyen d’une 
notation vectoridle simple. Designons par les lettres 


santes 


e,, 
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Les quatre derniers sont des vecteurs d'espaee. Avee ecs notations 
la « quantite de mouvenient-inasse » d'un point materiel de masse m 
est representee par 

/ dx dy dz \ 

Si nous convenons encore de designer par une lettre m un point 
d’Lnivers (/, .r, v, - la derivee — de ce point par rapport an temps 
est le vecteur dTnivers dc couiposantes 


i. 


dx 

dt’ 


dy 

dt ' 


dz 

dt’ 


on voit quo la « quantite de mouvement-masse * d’un point materiel 
est representee par la notation 

eh n 
m — — • 
dt 

Les points et les vecteurs (lihres) sont des formes geometriques du 
premier degre. On peut eonsiderer aussi des formes geometriques du 
second degre, qui represenlent des systemes de vecteurs glissants. 
On desigue par [mm j le vecteur glissant qui a pour origine le point 
d’Univers m et pour extremite le point dT nivers m'. Ce vecteur 
glissant a dix coordonnees pliickeriennes qui sont les determinants 
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K. C.AKTAX. 

du deuxienie ordre formes aver l«* Tableau 


t l i r z, 

i f' .c v' z : 

on a evidemment 

[m'm | = — [mm ]. 

De ineiue on desi;mera par jmej le verteur "lis>ant obtenu on por- 
lanl a partir du point d'l nivers m mi verteur «*<| u ijiol lent a on verteur 
doniie e; li*s di\ coordnnnees pluckeriennes do re verteur *;lis>;i tit 
sont formers avre 1«* Tableau 


• < ■>' y -• 

O '< i r, 

oil li^nmit dans la semiide 1 iirii*- les mmposantc- du verteur e. Enfin 
la notation |ee | desi"in*ra le hixTcteur dont les dix coordoiinees sont 
formers aver le Tableau 

o > ■ I, r. 


des roinposaules ties deux veetenrs libre> s et e . 

Dans ehaeun des cas precedents le verleur "li>-ant mi le biveetem* 
eonsidere pent el re rrj'arde rniiuiie le produit ( exterieur ) des deux 
farteurs. <| it i sont des formes "roinetriqurs du premier de^rr ( point 
nu verteur libre ). Le produit dr deux formes ifrometriques quelron- 
rjues du premier de«re salislait it la b*i distributive. mais change tie 
si^ne aver I’tirilre ties farteurs. 


12. Le verteur jilissanl t|iii a pour ori"ine le point tlTuivers m qui 
represent!* mi point materiel doom* a un instant donue, et qu'mi 
obtient eu portant a partir tie re point sa <> quantile do inoiiveiiient- 
masse », a pour expression 

I >/m I 
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oil F represente le \erteiir irlissant « force », eontient, i*n mime temps 
que le principe fondameiital tie la Dvuamiqur. le theoreme des 
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moments cinetiques: elle condense en efTet les dix equations 
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13. Revenons a la Mecaniquc des milieux continus. Designons 
par G le vecteur *rli>sant qm represente I'element a trois dimensions 
de « quantile de mouveincnf-masse » et par F le vecteur "lissant qui 
represente la force de volume elementaire. Les equations dr, la Meca- 
nique sont condcnsces Hans la forrnulc 

<<>> G'=[rfiF]. 

On a du reste ici 


G = [me u ]U - [me, 1 11,-+- [me.]^ - | me,]II.. 

F = | met ]\ dr dy dz - [me.]t i/j-i/n/: + [ae,]Z(//(/i </;. 

Le calcul de G' pent se faire en tenant cmnpte de (’equation 
dm — e, dt + e, d.r -+- e, dy — +- 6j dz: 

il donne 

G = lme,]ir+ f me, ] II j -+- [me t ]II', -+- [me, [III 

-+- [e 0 e,][v//ll, - 'Aril] -+-[«„«.] ['/Ml_. — '/. II ] - [e 0 e, j fo7II : — rf;!!] 
-t- f®t •»][*/>' IL— dz\l y \ -+- [«,e,] [f/dl, — 'Aril.] |e,e-][</.» II, —dt 11^]. 

4mm. Ec . Norm. % (1), XL. — NufUUK »y»!i. 44 
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Chapter III 
Attitude Estimation 

A new concept for estimating the three angles which specify the 
orientation in space of a rigid body is now presented. The estimation pro- 
cedure is analogous to phase -locked loop phase estimation: an observed 
function of the unknown angles is modulated by a function of the estimated 
angles, the resultant function is filtered by a linear time -invariant system, 
and the system outputs are the angle estimates. Two separate methods for 
implementing this concept are discussed, one based on direction cosines, 
the other on quaternions. No attempt is made to show that the particular 
estimation method presented here is an optimum method, in the sense of 
minimizing an error criteria. The aim of this chapter is to present a work- 
ing technique that offers potential savings in hardware, in certain cases, 
over, for example, methods based on Kalman filtering (see Jazwinski( 3 1) ). 

Before beginning the actual description of the attitude estimation 
procedure developed in this research program it is perhaps instructive to 
survey those features common to any estimation problem. Those features 
are three in number: 1) a parameter space, a space on which the parameters 

to be estimated are well defined, 2) an observation space, where functions of 
the parameters, corrupted by noise, can be observed, and 3) two mappings, 
one from the parameter space to the observation space which defines the 
functions to be observed in terms of the parameters and noise, and one from 
the observation space to the parameter space, which defines an estimation 
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procedure, a method for processing observations in order to estimate the 

J 

desired parameters. For the problems to be discussed here, the parameter 
space is a real three dimensional Euclidean vector space in which three 
angles take on their values. The observation space is a real Euclidean vector 
space. The mapping from the parameter space to the observation space can 
be done in many ways, but only two will be discussed here, a direction cosine 
mapping and a quaternion mapping; each of these mappings may be viewed as 
an exponential mapping of the three angles from the parameter space to the 
observation space, loosely speaking. For both mappings, observations are 
functions of sines and cosines of the unknown angles, as well as noise. The 
estimation procedure, the mapping from the observation space back to the 
parameter space, is the subject of this chapter. 

In order to assess the worth of the estimation procedure, one approach 
would be to define an error criteria, a measure of how much the actual 
parameters differ from their estimates. The error criteria considered in 
this chapter, often implicitly, is the difference between the actual value of 
an actual value of an angle and its estimate. An optimum estimation procedure 
minimizes the error criteria to its smallest possible value, assuming the 
error criteria actually has a minimum and that an optimum estimation rule 
actually exists; the question of optimum estimation will not be addressed in 
this chapter. 

3.1. Attitude Characterization 

This section discusses two methods for specifying the spatial orienta- 
tion of a rigid body, direction cosines and quaternions. The material presented 
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here is largely tutorial; Goldstein(lO) is an excellent reference for material 

i 

on direction cosines, while Whittaker (28) is fine for qurternions. 

3.1.1. Direction Cosines 

3 

Consider two right handed cartesian coordinate frames in II with 

the same origin labeled A and B. Assume that when frame A is rotated 

about some axis, eventually it will coincide with frame B (see Euler's 

Theorem -Whittaker, p.2). If an arbitrary vector has coordinates 

r A = (x,.y 4 ,zj in frame A. and r = (x , y_ , z_ ) in frame B, then it i 
— A AAA id id Jd x5 

well known 


is 


3.1) 


B 


B 


"B 


d ll d 12 d 13 


21 


d 22 d 23 


d 31 d 32 d 33 


r. = Dr , 
> — -A 


where D is the direction cosine matrix from frame A to frame B. The 
physical interpretation of the elements in D is simple; e.g., d ^ is the projec- 
tion of a unit vector along the z-axis in frame A onto the x-axis in frame B. 
The elements in D are constrained by 


3.2) 


T ' T 

D D = D D = I 


where D is the transpose of D. Since there are six constraints in this 
matrix equation. 


2 2 2 
d ll + d 12 + d 13 


d ll d 21 + d 12 d 22 +d 13 d 23 " ° 


2 2 2 

d 21 + d 22 + d 23 = 1 

2 2 2 

d 31 + d 32 + d 33 ' 1 


d Il d 31 + d 12 d 32 +d 13 d 33 = ° 


d 2t d 31 + d 22 d 32 +d 23 d 33 = ° 


1 


while ^ contains nine parameters total, D can be characterized by three 
independent parameters. These three parameters are angles, denoted by a 
three - tuple a_ = (a ,a ,a ), represent rotations about x,y, or z axes: 

3.3) >= EXPKfc,) t*P{ EXP( k=|,3,3 L,* 



0 0 0 


0 0 1 


0-10 

L = 
— x 

0 0-1 

0 1 0 

L = 

— y 

0 0 0 

-1 0 0 

L = 
— z 

1 0 0 
0 0 0 


* «■ 




_ w 


In all cases, D is specified by the three independent angles a ,a ,a . 

1 4 j 

However if either 



i) ki * i-j 

4*»D 

'r- 

V K 

fc= 0,±l,±J,.- 

or 

ii) t. 

hub 





then D can be shown to be composed of sines and cosines of (a^+a^) or sines 
and cosines of (a^-a^), depending on the exact choice of L and L in terms 
°^|± x , Ly- In this case D is "singular," in that specifying only two 

independent parameters, a^ and either (a^a^) or (a -a^) un iquely determines 
D. 


If coordinate frame A is rotating relative to frame B, then the direc- 
tion cosines from A to B change with time. The derivative of D with respect 
to time is 



I 4|Jr, + 


Q.i, . , 

«< * + 


ttib» 

v t *• A 
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where w is called the angular velocity of frame B with respect to frame A. 
3.1.2. Quaternions 

A second way of characterizing a rotation in space is by a quaternion. 
A quaternion q may be viewed as a four-tuple (with one constraint) 

3.5) + + d(l) 

A A A 

where i, j, k and 1 are unit quaternions, which multiply according to the 
rules of quaternion multiplication: 


II 

A A A 

-j i = k 

A A A 

1 i = i 1 =i 

a2 

i 

*2 

= J 

!V> 

II 

A A A 

-k j = i 

A A A 

1 j = j 1 = j 

i 2 

= 1 

H-> 

II 

1 

II 

A A A 

1 k = k 1 = k 




The conjugate of a quaternion q, denoted q , is defined to be 

3.6) q = -ai -bj -ck + d(l) 

The four parameters of a quaternion are constrained: 

3.7) qq = qq = a +b + c +d =1 

Since any rotation in space is a rotation about some axis, a quaternion may 
be viewed with respect to a reference coordinate frame as 

3.8) C60(?) + £*(!?)[ K ! j 

where at any instant of time the rotation is about a unit length vector, called 
the instantaneous axis of rotation, and 



COS A 


projection of instantaneous axis of rotation 


onto x-axis of reference frame 
cos^ = projection of instantaneous axis of rotation 
onto y-axis of reference frame 
cosY = projection of instantaneous axis of rotation 
onto z-axis of reference frame 
w = amount of rotation (in radians) about instan- 
taneous axis of rotation 

Using the power series definition of an exponential plus the rules of quaternion 
multiplication, this can be rewritten as 

3.9) ^ = EXp{ jL + J 

If w = w +2nn, physically this represents a rotation of w radians about the 
o o 

instantaneous axis of rotation. Note that 

^ (§)C t <***+ ] 

3.10) (-0* [ ttf( **) + j fc 1 j - (-0 

that is, both (if n is even) and -q^ (if n is odd) represent the same rotation 

in space. This fact is reflected in how quaternions transform vectors in 

frame A to vectors in frame B. If an arbitrary vector has coordinates 

r . = (x . ,y . , z . ) in frame A, and r = (x , y_ , z_.) in frame B, where 
— A AAA — ±4 d d a 

frame A and B are related by a rotation about some axis, then 
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Since this can also be written as 
3 . 12 ) r B = (-g) r A (-g + ) 

a quaternion is a point on a sphere in Euclidean four- space, but the rotation 
that quaternion represents must be considered as two antipodal points 
(a, b, c , d) and (-a, - b, - c , -d)^ on that sphere . 

If coordinate frame A rotates relative to frame B, then the quaternion 
relating those two frames changes with time; it can be shown (Whittaker (28) , 


p. 16) 


3.13) 


dt 


J. 
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z u X 

0 -w w 
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where w is the angular velocity of B with respect to A. 

3.2. Attitude Estimation With Direction Cosines 
3.2.1. Deterministic Attitude Estimation 
3.2. 1.1. Model 

A block diagram for a system which estimates the three angles which 
specify the orientation of a rigid body from direction cosine measurements 
is shown below (without any noise sources): 
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It is assumed: 

1) all components of D^, the received direction cosine 
matrix, are observable 

A 

2) and D, the received and estimated direction 

cosine matrices, respectively, are each described 

by three roll -pitch -yaw angles, a and a, respec- 

R 

tively 

3.U) D r = exp (x r L^) exp (y R I^) e*p { z r LJ a R = (x R , y R , z R ) 

3.15) D = exp (xL ) exp (yL ) exp (zL ) a = (x,y, z) 

—x y z 

where 



0 

0 

0 


' 0 

0 

1 


~0 

-1 

o~ 

L = 

— X 

0 

0 

-1 


0 

0 

0 

> -z = 

1 

0 

0 


0 

1 
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-1 

0 
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0 

0 
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A 'p 

results from postmultiplying D R by ^ , the transpose (denoted with 
a "T'Jof D. . 


3.16) E = D r D T = 


E 11 E 12 E 1 3 

E 21 E 22 E 23 

E 31 E 32 E 33 


If each of the six angles is small compared to one, then E can be approximated: 


1 



U p l«l 
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The 3x1 error 

vector g may be defined as 


S s 

E Si <*,) - £» Cl-M,) 


3.18) *2 ' 

E » t«|l - 

0 l 

V' 


0 I 


Each error component may be chosen in many ways; it is not obvious that 
one choice is better than any other choice. Only two will be considered 
here, to provide a concrete example for discussion, and to simplify algebraic 
manipulation: 


S = 


«.* - 


3.20) " E *‘ 

The error vector e^ is used as*an input to a linear time -invariant system which 


3.19) *3 ’ ^*3 

* £-> 


generates angle estimates, a. 


3.21 


£ = A a + B e 

A „ 

a = C a 


a (t ) = a 
— o — o 


where A,_B, and are 3x3 matrices, and a_ are 3x1 vectors . For the two 

examples to be considered, A = so this system is a pure integrator. 

T 

The angle estimates are then used to generate E) , closing the loop 

T 

The intuitive operation of the loop is simple: postmultiplying D by D 

R 

undoes the rotation represented by D ; the negative feedback of the loop 

R 

eventually zeroes £, eventually making E the identity matrix and a = a 

R 

For future reference the off diagonal elements of E are now explicitly spelled 


out: 
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3.22a) clx) &lfl) -Sl^r) Steely) 

3.22b) Eijc -tl«JrH(^r-|)cc{)slJ)- ft^r) +• "tyrW) tO$) 

3.22c) E Jf : ttj) jsCXr)3t3rlcf3r'^ +Ct3J r )s(^-3')j -S(KrH(^ r )5Cj) 

3 . 22d) Ejj = - tt<) «t j3 { il^r'lsUr) t(} r -|) + ClxJ^r'}* J " 5C *»0 r) cCjc) tty) 

+ «X){ c^ r H^ri)“*tXr)^ r U(2 r ^} 

3.22e) E it = c(^{«X r )*^ r -j)-*t^ec*r)etV'5^ * eUr) ety) « j) 

3 • 22f) E^ = itxlity) |stx r )s(j r -^ -Sty^) ctvOctjr'^l " 5£ x)tly) cU r ><ij r ) 

<-C(x) [ 3(»jr) ttXr) tC 3r'3^ 

where sine and cosine have been abreviated to s and c, respectively. To 
carry the analysis further, two examples are now examined in detail. 

3.2. 1.2. Example 1 --Constant Orientation in Space 
Throughout this section a^ is assumed constant but unknown. If 



o o o" 


1 

x°“ 

o 

o 

1 


1 o o~ 

A = 

0 0 0 

B = 

0^0 

c = 

0 1 0 


0 0 0 


1 

o 

o 

09 

r 


0 0 1 


3.23) a = B«e (a_,a) a(t ) = a 

— — R — — o — o 

If (3.19) is used for e, these three equations can be written out 
explicitly: 

* ' &x{ S <*) Sl^[s(^ 5 (^ r - \) -s(^) c (X/)c(} r - 3 l ] - s(i)®^ 6(X^c(^) 
♦ o(x)£s(^c(x^sl}r-}) f SCXrHC^-^)) \ 

5~ ^l-clj r ) c ( v -^cCx)s(5^- c ^5(^-J)sCx) + styOcCxWpJ 

V S\ ct ^H^s c ^o(y^ + c(kr‘) 4 (^')]-su r ) 6 (^s (^5 


3.24) 
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The figure below shows a sampled-input, discrete-time approximation to 
these equations 
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The figures on the following page show the results of a computer simulation 

of the equations above using the discrete time filter (Figure 3.3 

a = 0, a = (60°, 45°, 30°); Figure 3.4 £ = 0, a „ = (-60°, -45°, -30°) ) 

— o — R — o — — R 

with A t = 0.1 seconds in each case , and B = B = B = 1 . 

x y z 

The next step is a linearized analysis of these equations. If a_ = 
(60°, 45°, 30°), and if it is assumed 

X ■= 60 ° * Sx 1 \« i 

3.25) ^ = 4 5°+ty U})«1 


V" 50%^ IJjki 


then using the approximations sin a * a, cos a * 1 for , and keeping 

terms to first order only, it can be shown that 



FIGURE 3. 3-3.4 

DIRECTION COSINE ATTITUDE ESTIMATION- -COMPUTER SIMULATION 
ESTIMATED ANGLES VS TIME 
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S% = -*K 

3.26) 

4 i = -h 


Note that estimation errors can be nulled out faster than in Figure 3.4 (where 


B = B = B =1) by choosing B 1 , B » 1, B » 1. A word of caution is 
x y z x y z 

in order: strictly speaking, this linearized analysis is valid only for 

while Figure 3.4 shows that near t = 0 there approximations 
are not valid. As long as estimation errors are less than 10°, i.e., t> 6 
seconds in Figure 3.4, the linearized analysis provides good agreement with 
simulation re suits . 

If instead of equation (3. 19), equation (3.20) is chosen for e^ , the 
results of both the computer simulations and linearized analysis are identical 
to those just discussed. 

The stable states of equation (3.23) were investigated next. 

First, the steady state solutions to (3.23) were found; if (3.19) is 


used for e^, the steady state solutions are 

$ *-X r + axt 

* , b£\)6M 


X*. > <\TT 

* b ExJfo) 


3.271 ^ 




C ODD 


X 13 v+au 

a boiD 

o odd 


Jr+c.it V* t7r J* 

where (a,b,c) are integers, while if (3.17) is used for e_, the steady state 


solutions are 


X - X r 4 ft* 


x e 


X a- -x r + air 


3.28) ~ bT 


bew £ =. .^Im, b ^ J B V + b * 


C.BMEM 


a odd 
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As a check on these solutions, a was held fixed, and the M.I.T. 

.K. 

Computation Center Fortran IV Subroutine "NONLIN" was used to solve 
e (a , a) = 0, with a grid of initial trial solutions for la. No more solutions 
were found in this way than those above; clearly, this is no guarantee all the 
steady state solutions are known. 

To illustrate this point, a different choice of e was made: 

«- x = 

3.29) 

A paper -and- pencil solution to £ = _0 yielded results similar to those in 

equations (3.27) and (3.28). However, using this <e , when a = 60°, 

a = 45°, a = 30°, "NONLIN" showed that e = 0 when a = 188.88°, 
ry rz — — x 

a = 209.75°, a = 39.13°, as well as at the expected solutions. 

y 2 

The stable states can be found by linearizing the equations of motion 
about the steady state solutions above, and then perturbing the equations 
away from the steady state solutions to see if solutions to the linearized 
equations are stable or unstable. To be more explicit, denote a steady state 



order terms: 
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3.31) 
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V*»S»V 3 h ,x M‘i«»v^ ♦ + 



Then in terms of perturbed variables (&()£%£%)* equation (3 . 10) becomes 



i 
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k *• 
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3.32) i 
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aj*l_ 
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where e (x ,y , z ) = 0 and "ss" means "evaluated at x = x , y = y , 

— ss ss ss — ss ss 

z = z ." Note that the eigenvalues of the linearized system can be deter- 
s s 

mined by the steady state solutions. For example, if (3.19) is used for e^ 


one set of steady state solutions is 

4 » x»- + ** 

^ = ^TT b Obbj C. Ob0 

v v- * 

It can be shown that ifx = y =z = 0, none of these steady state solutions 

r r r 

are stable, but ifx = If , y - z =0, then a stable steady state solution 

r r r 

exists for a,b,c all odd. The figure on the next page shows the result of a 


computer simulation of equation (3.24) with = 1, and 

X* » 0° *r> = 

V°° V *> # 

Clearly the stable steady- state solution to this simulation is 


a = -i=o° = 

Y r + (-1) 7t 

$ = is* = 

-‘Jv * 00 ir 

r- -vr - 

\r + 


Degrees 
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DIRECTION COSINE ATTITUDE E S TIMA TION 
COMPUTER SIMULATION RESULTS 
ANGLE ESTIMATES VS TIME 



FIGURE 3. 5 
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3.2. 1.3. Constant Rate of Change of Roll- Pitch- Yaw Angles 


Throughout this section, a and £ are assumed -o be of the form 

t4< r » J* ■'»**!* ] a r = it*i r 


* r J = ifjt+ 


V 


& = i*+t 


where f = (f ,f , f ) are known and constant, *P =(^ , ) are unknown 

— x v z T rx ry rz 

^ a 4 ^ 

and constant, and *{=(«$ ,<P ,<€ ) are estimates of = (*£ , K? , ^ ). The 

x y z r rx ry rz 

reader is cautioned not to confuse this set of circumstances with a constant 
angular velocity; from Goldstein( 10) , angular velocity in terms of roll-pitch- 


yaw angle s can be shown to be 


0 cos a -sin a 


0 sin a 


0 cos a - sin a 

rx rx 

0 sin a cos a 


cos a 


0 sin a 


1 0 


cos a -sin a 0 cos a 


a + 

r y 


Only under special circumstances (e.g., f^ = f = 0) does constant rate of 
change of roll -pitch -yaw angles reduce to constant angular velocity. 


The system block diagram is slightly modified to take this into account: 
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The only change from the earlier system is a processor which demodulates 
and filters • ,fc) to give . If £ is given by (3 . 16) , then 

a possible processor is: 
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Each low pass filter is assumed ideal, passing signals at zero frequency 
without attenuation, but perfectly blocking signals outside this band from 
passing . 


The demodulated and low pass filtered error, e/ , is used as the input 
to the linear time invariant system: 

< = A 3 +• 6 e' 5 ( 0 = ? 0 

3.33) A ^ 

= C(f 


For the example above, e’ is 


4* = g + s ( dj, + + iJsUa+cJ-jl 

£ dy) A 


. * A. 

<«= 8 


4^ = $ L$( A- S( 






The linearized analysis of these equations is straightforward: if *£ = 0, 

A r A 

then d = and if | cc l j Kjl^l , **■ I i then to first order in 
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3.34) 
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a 

4 } , t 11 *'”” 


t*t 0 


That is, for small initial errors this system will estimate the correct phases. 
The steady state solutions can be discussed for this example just as in the 
previous example. No computer simulation of the equations of motion was 
carried out. 


3.2.2. Stochastic Attitude Estimation 
3.2.2. 1. Model 

A block diagram which includes noise sources is shown below, virtually 
identical with that in Figure 3. 1: 
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Two noise sources are evident: 

1) a , noise generated by sensor measurements, 

— n 

assumed added to the true angles a 

— s 

2) n, noise generated by modulation and amplifica- 
tion of sensor measurements. 

The received direction cosines are: 

3.35) D = exp [ (s + n' ) L ] exp [ (s + n 1 ) L ] exp [ (s + n' ) L ] 

— R x x — x y y — y z z — z 

i.e., products of exponentials of random processes. Ito (36) and McKean (41) 

have done some preliminary work in characterizing this particular process. 

However, the problem considered here , finding the aposteriori probability 

density of a, is apparently still unsolved. From here on, a is assumed 0, 

— — n — 

while n is assumed to the only source of uncertainty in estimating a given 

~~ S 

D . This is an ad hoc assumption, done only in the hope of making the new 
K, 

problem analytically tractable. 

3. 2. 2. 2. Fokker-Planck Analysis 

For example one, including amplifier noise n but not sensor noise a , 

(3. 23) becomes 

* _ 

x = B e (a, a = a ) + n 
xx r — s x 

* T 1 

y = B e (a, a = a ) + n E(n) = 0,E[n(t)n (s)] = — £ (t-s) 

yy r— s y — 2 

• 

z=B e (a, a = a ) + n 

z z r — s z i 

where n is assumed to a zero-mean white Gauss-Markov random process, 
uncorrelated between axes. 


N 0 0 

x 

0 N 0 
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0 0 N 
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It is straightforward to derive the Fokker-Planck equation for the 
aposteriori density function of S(t), given <5 + n from an Initial time t to the 
present time t (see Viterbi (48)): 


This equation has an initial condition, in that a(t ) is known: p(a(t ) ) = i (a - a ). 
Since e is periodic in each component, p has a boundary condition: 

p(a(t) 1 n, [t , t) ) = p(a(t) + (a 2tt , b2w , c2f )J e, n, [t^.t) ) } a,b, c integer s 


This boundary condition suggests a possible approach to either solving 

the Fokker-Planck equation or approximating a solution is to expand the 

probability density in a Fourie’r series for each of the three estimated angles; 

this approach was never seriously investigated. 

Either equation (3. 19) or (3.20) can be used for e ; in neither case was 

it possible to solve the Fokker Planck equation explicitly. 

For example two, again including amplifier noise n but ignoring sensor 

noise a , (3.33) becomes (B~l) 

— n 

• ' 

$ = 2s(d -d -d )-2s(d -d +d )+s(d +d )+2s(d +d )+2s(d -d )+n 

# x yxz yxz x y xz xzx 


3 . 38 ) 
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= -s(d +d )-s(d -d )+n 
y y z y z y 

= s(d +d )+s(d -d )+n 
z z x z x z 


E(n) = £ E(n(t)n T (s) ) = ~ & (t-s) 


N 0 

x 

0 N 

y 

o o 


o 

o 

N 


63. 


The Fokker-Planck equation is identical in form to that in equation 

! 

(3.37), with identical boundary conditions and initial cciditions. Again, it 
was impossible to solve this equation explicitly for the aposteriori probability 
density function of a. 

3.2.2. 3. Angle Estimate Skipping 
Recall that without noise present many steady- state solutions exist 
for the two examples discussed. With noise present, it is possible for the 
systems discussed here to skip from one steady-state solution to another, 
with the skip caused by a burst of noise; computer simulations were carried 
out to observe this skipping, and showed it indeed does occur. 

There are two ways for the estimated angles to skip: first, any esti- 
mated angle may increase by a multiple of 2n radians; second, all the angle 
estimates may hop from one type of steady -state solution to an entirely dif- 
ferent one (cf. the three different types of steady-state solutions in equations 

(3.27) - (3.28)). The first case corresponds to a rotation in space about 
« 

some axis by a multiple of 2tr radians; physically, this corresponds to no 
change in the estimate of orientation in space. The second case does corre- 
spond to a change in the estimate of spatial orientation. The only approach 
known for investigating these two cases is to check the stability of steady- 
state solutions in the absence of noise, as was discussed in sections 3.2. 1.2-3. 
Bounds on how frequently an angle estimate might skip are not known at the 


present time . 
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3.3. Attitude Estimation With Quaternions 

3.3.1. Deterministic Attitude Estimation 

Because of algebraic complexity encountered in using direction cosines, 
as well as the inability to solve the Fokker-Planck equation for the aposteriori 
statistics of a, the same attitude estimation method was implemented with 
quaternions rather than direction cosines. This implementation is now 
discussed. 

3.3. 1.1. Model 

A block diagram for a system which estimates three roll -pitch -yaw 
angles which specify the orientation of a rigid body from quaternion measure- 
ments is shown below (without any noise sources): 



MErMOft.'iL&S.S 


FUcoOfc 33 
oPlWAI'OfS SPAtfe 


It is assumed: 

1) all components of q , the received quaternion, are 

r 

observable 

2) q and q, the received and estimated quaternions, re- 

_ir 

spectively, are each described by three roll-pitch-yaw 

A 

angles, a^and a, respectively. 
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3.39) a r = [c(x r /2)+fs(y r /2)] [c(y r /2)+ j s(y t /2)][c(z r /2)+k s(z r /2)] 

3.40) £ = [c(x/2)+ i s(x/2)] [c(y/2)+ j s(y/2)] [c(z/2)+! s(z/2)] 

A A. A 

where (l,i,j,k) are unit quaternions which multiply according to the rules of 
quaternion multiplication. £ is postmultiplied by , the quaternion con- 
jugate of q (denoted with a "+"), which produces an error quaternion, q 

-*e 

• A A A 

3.41) q = q q r = e i+e^ j + e k + e 

e a r x ” z ^ 


w 


If each of the estimated angles as well as the actual angles has magnitude 

much less than unity, £ can be approximated 

3.42) e s (x - x)/2 e 3? (y - y)/2 e 2f (z - z)/2 e » 1 
xr Jy r ?zr w 

Unlike the direction cosine case, £ can be chosen here in only one way, 

e = (e , e , e ) . 

— x y z 

e = (e ,e ,e ) is used as an input to a linear time - invariant system which 

— x y z 


generates a: 


3.43) 


a = A a + B e 


1 = Ca 


A . 

a is then used to generate £ , closing the loop. The intuitive operation of a 
. * * 

quaternion phase -locked-loop is identical to that of a direction cosine 

• * . 

phase-locked-loop: postmultiplying £ ^ by £ undoes the rotation represented 
b Y S, ; the loop negative feedback action nulls £ to zero, eventually making 
£^ = 1 and = £. For future reference the components of £ are explicitly 
stated: 


f 
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3.44a) 

3.44b) 

3.44c) 


e = s(d /2)c(d /2)c(d /2)+c(d /2)s(s /2)s(d /2) 
xx y z x y z 

. * 

s = X + X 
x r 

, A 

d = x -x 
x r 

e = c ( s /2)s(d /2)c(d /2)-s(d /2)c(s /2)s(d /2) 
y x y z x y z 

, A 

s = y + y 
y r 

. * 

d = y -y 

y r 

e = c(s /2)c(s /2)s(d /2)+s(s /2)s(d /2)c(d /2) 
z x y z x y z 

A 

s = z + z 
z r 

d = z -z 
z r 


To proceed further, the two examples discussed in the direction cosine case 


are re-examined. 


3. 3.1.2. Example 1- -Constant Orientation in Space 
Throughout this section, a is assumed constant but unknown. If 
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then equations (3.43) and (3.44) can be combined: 

x = {s(d / 2)c(d /2)c(d /Z)+c(d /2)s(s /2)s(d /2)?B x(t ) = x 
*■ x y z x y z -> x oo 

3.45) y = fc(s / 2) s(d /2)c(d /2)-s(d /2)c(s /2)s(d /2)\B y(t ) = y 
*• x y z x y z > y o o 

• 

z = fc(s / 2)c(s / 2)s(d /2)+s(s /2)s(d /2)c(d /2)7 b z(t ) = z 
L x y z x y z i z oo 

This discrete -time sampled-input approximation in Figure 3. 2 was 

used to approximate these equations; the figures on the next page are the 

results of two computer simulations (Figure 3. 10, x = 60°, y = 45°, 

z = 30°, a =0° and Figure 3. 11, x = 60°, y = 45°, z = -30°, a = 0°) 
r — o — r r r — o — 

with At =0.2 seconds and B = B = B =1, 

x y z 

The linearized analysis, assuming a = (60° , 45° , 30°) , and «, = 

(x r +i x,y r + Sy,Z r +£ z) with first order in shows 


Degrees 


t 


QUATERNION ATTITUDI 


COMPUTER simulat: 



Time (seconds ) 


FIGURE 3.10 


Degrees 





FIGURE 3.11 
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S* s 

6, [-iix-fij] 

<hfe>) = 

3.46) S 


Sy It) z 

H * 


a 


-iVfc.-f) , 

t 




The steady state solutions to e = (e , e ,e ) = 0 are 

— x y z — 


x = x + n a 
r 


x = x + a IT 
r 


y = y + IT b a,b,c even; y = -y + bjf 
r r 


a,b,c odd 


z = z^ + Tfc 


z = z + C7T 
r 


Again, an analysis similar to that in Section (3 .2 .1 . 2) must be carried 

out to determine which states are stable and which are not. For example, a 

computer simulation of these equations was carried out for 

.o 


II 

0° 

X = 

120 

o 

A 

y = 

0° 

II 

> 

105 

o 

A 

„o 

r 

o 

z = 
o 

0 

z = 

r 

90 


and the results were similar to those in Figure 3.6, with the steady state 
(denoted ss) angle estimates being 


£ = -60' 
ss 

A •7C t 

y =75 
ss 

A c 

z = -90 
ss 


x r + ax 

-y r + M 


z r + cjT 


a,b,c odd 


3.3. 1.3. Example 2 --Constant Rate of Change of 
Roll- Pitch -Yaw Angles 

In this section a. = H + Sr, and a. = £_ t + if , where f_ is known and 

A 

constant, is unknown and constant, and ^ is estimate of . The block 

r “ r 


diagram of the attitude estimation system becomes 
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^ ®,' - ^ e x’ e y’ e z^ * s ^ n P u ^ to the linear time invariant system 



then the equations for are 

% = *X l j % 

3.47) ^ r ^ ^a)= 4 jj & 

"*» ~ ®ii c -W n ,^»)/a]e[[<« J i'< 9 )/j]s£t>< r} - 3 j)i!.]} lt»l c <? }tl 


The linearized analysis of these equations (assuming tf = 0, and l^l^l 

^ i * 

1^|^<M) is identical to that in (3.36). The steady state solutions, which 


are also the stable solutions, are 
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i r <l rJ i + 

X 

iu - 'tr* 4 ait b 

3.48) X TJ * 

^ ^ + *tt 


4,lo,G faST£<*6(3J 


3.3.2. Stochastic Attitude Estimation 


3. 3.2.1. Model 

A block diagram which includes noise sources is shown below. 


»\ 
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where a arises from sensor measurements, n from amplification and modu 
— n ~ 

lation as in the direction cosine case. 

The received quaternion is: 

A -4 

q = [c(s + n'/2) + i s(s +n'/2)] [c(s + n'/2) + j s(s + n' / 2)] 

-*r L x x x x y y J y y 

A 

[c(s + n' /2) + k s(s + n / 2) ] 
z z z z 

Again, no methods are known for handling a ; from this point on it will be 
° — n 
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assumed a =0, while n is a zero -mean white Gaussian random process, 
— n — — 

uncorrelated from one component to the next. 

3. 3. 2.2. Fokker-Planck Analysis 


For example one, including n in equation (3.35) shows 

t 

x = B £s(d /2) c(d /2) c(d /2) + c(d /2) s(s /2) s(d /2)"l + n 
xt x y z x y z J x 

* 

3.49) y = B fc(s / 2) s(d /2) c(d /2) -s(d /2) c(s /2) s(d / 2 )} + n 
V * x y z x y z ' y 

z = B fc(s / 2) c(s / 2 ) s(d /2) + s(s /2) s(d /2) c(d /2)}+ n 
zlx V z x y ZJ2 


where 


E [n(t) n T (s)] = t ! (t-s) 


N 


x 

0 


N 


y 

o 


N 


1 E[n(t)] = 0 


The Fokker-Planck equation is straightforward to derive: 

9 • 0 rn . f n i Mi 


3 . 50 , 

_ ir s i£ * ^ ^ £n 

*jL» jS** * j Tf-I 

The initial condition and boundary conditions are 






Three assumptions are made 


i) a = a =0, the actual angles a are all zero 
— r — s — — r 

ii) N = N = N = N , the noise covariance is identical 
x y z 

in all three axes, 


iii) B = B = B = B, identical gain in all three error 
x y z 

amplifiers 
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then the steady- state Fokker- Planck equation can be written quite simply: 

3 . 52 , ^ 

© © <D 

If the solution is assumed to be p^, where 

3.53) Pj = C exp -|r [c(x/2) c(y/2) c(z/2) - s(x/2) s(y/2) s(z/2)] 

C = normalization constant 


then the first and third terms in (3. 52) are zero, but not the second. On the 
other hand, if the solution is assumed to be p^, where 

3.54) p 2 = C exp ^ ^[c(x/2) c(y/2) c(z/2) + s{x/2) s(y/2) s(z/2)] | 

C = normalization 
constant 

then the second term in (3.52) is zero, but not the first or third. Even 
though the solution (3.52) seems within reach, no solution was found. 

For example two, including amplifier noise n, (3.47) now becomes 

i 3 «, j c J « v V* +,, > 

3.55) . + V *1*i**> 

^ 8j{c(^}eC 
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The Fokker-Planck equation is identical in form to (3.50). Again, assume 

tf = 0, N = N = N = N, and B = B = B = B, so that the steady state 
r — x y z x y z 

Fokker-Planck equation becomes: 

3 - 56 » k** < itifoutytiipn * ^ + 

p- pell 


The solution to this equation is 

3 - 57) f til V) - C-b*p| bl'L'ijb)! 


where C is a normalization constant, chosen such that 

r ^ 

l \ 


, a* 

j/ 


3 - 581 \ \ J. pCil* ri ft.'<t-»l3 A«MjA 

V-^ir w 


- 1 


if 


'fy" 1 * 


p S C. EXP J ^ 


a a a j 

A A 


3.59) 


which appears Gaussian with variance (N/2 8^. 

To keep the steady- state error variance small, B should be chosen such that 

Recall that in the linearized analysis of the model, B corresponds to bandwidth. 
Speaking intuitively, the larger the bandwidth B the smaller the steady- state 
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error variance; this is analogous to the situation in an angle modulation com- 
munications system, where by increasing the system bandwidth noise effects 
are suppressed. 

At this point, the steady-state Fokker-Planck equation has been ex- 
plicitly solved for a special set of circumstances, a quaternion implementa- 
tion of a particular estimation procedure. Why not return to the direction 
cosine implementation of this estimation process and solve the steady- state 
Fokker-Planck equation under the same set of circumstances as was just done? 

In fact, this was attempted, but the attempt was unsuccessful . The 
reason for this failure is apparently fundamental, based on the observation 
that the observation spaces for direction cosines and quaternions are essen- 
tially different vector spaces. Recall the underlying parameter space for both 
estimation procedures is a real Euclidean three-dimensional vector space. 
However, the direction cosine observation space and the quaternion observa- 
tion space are real Euclidean vector spaces, of dimension nine and four, re- 
spectively. Next, elements in the observed direction cosine matrix and 
components of the observed quaternion are functions of sines and cosines of 
the three angles to be estimated; the explicit relationship between each observed 
quaternion component and the elements in the observed direction cosine matrix, 
or between each element in the observed direction cosine matrix and the 
observed quaternion components is straightforward to work out but algebraically 
complex (this relationship reflects the fact that, loosely speaking, a quaternion 
when squared become s a direction cosine). Finally, since the observation 
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spaces are different in each case, the actual details of the estimation proce- 
dure are different; for the particular example considered here, the demodu- 
lators differ for the direction cosine case from those in the quaternion case. 

If it could be shown that each estimation procedure was an optimum procedure, 
then perhaps it could be shown that the two estimation procedures were 
equivalent; however, the rule for estimating angles discus sed here is a 
heuristic one, and it is not at all clear what the optimum estimation procedure 
should be, so it is not surprising that the two demodulators are different. 

In summary, while it was fortuitous that the Fokker- Planck equation 
could be solved under a special set of circumstances in the quaternion case, 
it is not at all clear why this good fortune should carry over under the same 
special set of circumstances to the direction cosine case. 

3. 3.2. 3. Angle Estimate Skipping 

The next page shows the results of a computer simulation of example 
one with amplifier noise included (equation 3.49), with a = (60°, 0°, -60°) . 

Note that two of the estimates of these angles are initially close to the correct 
estimate, and then skip by 2 71 radians, while the third estimate always remains 
close to the correct angle. Some of the qualitative issues as sociated with 
estimates of angles skipping have been discussed earlier (3. 2. 2. 3), and will 
not be repeated here. No bounds on how frequently the estimates might skip 
are known at the present time. 



Angle estimates (degrees) 


Figure 3. 14 


An Example of Angle Estimate Skipping --Constant 
Orientation in Space- -Quaternion Implementation 
of Attitude Estimation Method 
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3.4 . Summary 

This chapter developed a new method for estimating the three angles 
which specify the spatial orientation of a rigid body, given noisy sensor 
measurements. Two ways to implement this method, one based on direction 
cosines, the other on quaternions, were discussed for two specific examples, 
fixed orientation in space and constant rate of change of the three angles. The 
theoretical deterministic performance limits of each implementation were 
covered: first, an analysis of the dynamics of each implementation (assuming 

samll estimation errors) was carried out, and second, the steady state and 
stable equilibrium points of each method were discussed. When noise was 
included in each implementation, the theoretical stochastic performance limits 
were quite difficult to pin down; in particular, while Fokker-Flanck equations 
could be derived for angle estimation error a posteriori probability density 
functions under certain simplifying assumptions about the nature of the observa- 
tion noise, these equations could rarely be solved. However, for a special 
set of simplifying assumptions, with a quaternion based method for estimating 
the unknown phases of the three angles which change at a constant known rate 
with time, the Fokker-Planck equation was explicitly solved for the steady 
state error probability density function. 
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Chapte r F our 


Solar Pressure Attitude Control 

4.1. Introduction 

A brief overall description of two possible synchronous orbit communi- 
cation satellites is now presented, before proceeding to a detailed description 
of the attitude control dynamics of each design. Both satellites orient them- 
selves in space using torques generated from solar or sunlight pressure. 

The sketches on the next two pages show two possible designs for an 
air traffic control satellite that would be in synchronous orbit over the North 
Atlantic. Each design has two sets of solar cell panels and a central 
antenna-electronics body, with three solar sails attached to each solar panel 
for attitude control. In both designs the antenna rotates about its central 
shaft every twenty four hours, minus a factor due to the earth's annual motion 
about the sun; in Design I the solar panels rotate about their shafts once each 
year, while in Design II the panels wobble from plus 23.5 degrees (with 
respect to the antenna shaft) to minus 23.5 degrees and back again once each 
year . The attitude control problem is to point the antenna earthward and the 
solar panels sunward, aligning two prescribed body axes which is equivalent 
to three axis control of the solar panels and a prescribed motion of the 
antenna with respect to the panels. 

The technological constraints involved in these designs are now 
sketched extremely quickly. First, thermal design is simplified for the 
solar cell power plant since nominally only one surface faces the sun. 

Second, Design I involves a set of bearings for each panel which must rotate 
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FIGURE 4. 1 
DESIGN I 


Nominally perpendicular 
to the ecliptic 
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FIGURE 4. 2 
DESIGN II 
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NOTE ; Solar panel 1,2; Rotate ^ 
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oscillation per year. 1 

Antenna electronics : One 1 

complete revolution per day. J 
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a axis of rotation 
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Solar panel 1 
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extremely slowly in space while passing DC electric power from the solar 
cells to the antenna electronics; Design II circumvents this potential problem, 
because its solar panel need only be a spring capable of plus -minus 23.5 
degree motion. Third, the power systems for both designs contain solar 
cells and batteries (for use when the satellite is in earth shadow- - seventy 
two minutes is the longest shadow duration, twice a year (Goddard (56)). 
Fourth, structural problems are complicated by possible sail flexing and 
bending. For a much more thorough coverage of these and other techno- 
logical problems, the reader is referred to the bibliography (Abel (53), 
Fleischer (55), Goddard (57), Likins (58) MacLellan (59>, M.I.T. (60)). 

This report will concentrate on the attitude control dynamics of Designs I 
and II from here on, not on the systems engineering aspects of communica- 
tions satellites. 

4.2. Derivation of Equations of Motion 

Assume the solar panels and antenna (in Designs I and II) are ideal 
rigid bodies with massless interconnecting shafts and solar sails, and the 
solar panels are identical and rotate in unison about their respective shafts. 
Fix an inertial right-handed coordinate frame at the center of the sun, with 
its z-axis normal to the ecliptic and its x-axis pointing at the center of mass 
of the earth once a year, at the winter solstice. From the previous discus- 
sion on Newtonian mechanics, the dynamical equations of motion for the 
center of mass of either design Are 
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4 -D f t s - s ft M=0 

4.2) 4- (tP -MR ) = tF 

dt — — cm 

where P= momentum of center of mass in inertial coordinates 
M = total satellite mass 


F = external forces acting on center of mass in inertial coordinates 
The center of mass kinematical equations of motion for either design are 


4.3) 


d_ 

dt 


R 

— cm 


where R 

— cm 


= V 
— cm 

position of center of mass in inertial coordinates 


V = velocity of center of mass in inertial coordinates 
— cm 

Quite an extensive body of literature exists for these equations, in classical 
celestial mechanics as well as modern guidance and control theory. It was 
felt any work done here would not be a significant contribution to this area. 
From this point on, solar pressure attitude dynamics will be the main con- 
cern, with center of mass dynamics left for future work. 

The angular momentum equations of motion for either solar cell 
panel are: 


— (I w) + wx(Iw) = N +N 

dt — p — p — p — p — p —sails — ant -►panel 


where I 

-p 


= solar panel inertia tensor in body-fixed coordinates 


w^ = panel angular velocity with respect to inertial coordinates 


N 


sails = sa -il torque on panel due to attached sails 


N , = antenna reaction torque on panel 

— ant -►panel 
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Similarly, the angular momentum equations of motion for the antenna are: 

4.5) ^-[l (w +w)] + wx[l (w +w)]=N 

dt — ab — ab — p — p — ab — ab — p —panels -►ant 

where I , = antenna inertia tensor in body-fixed coordinates 
— ab 

w = antenna angular velocity in body-fixed coordinates 
-“ab 

N , = total panel reaction torque 

— panels -► ant 

Adding these two equations, and using the fact reaction torques are equal and 
opposite, the total angular momentum equation of motion is found to be 


4.6) §T [d ^ + 21 )w + I 


w 


] + w x[ (I + 21 ) W + I w ] = N 


dt — ab — p ~p — ab — ab -~p — ab — p — p — ab — ab J — sails 

where N , = total sail torque generated by all six sails 

s ail s 

A final kinematical equation is needed to describe how the direction cosines 
from body-fixed coordinates to inertial coordinates evolve with time: 


, d _b _b ,, T b 

4.7) — D. = D. W. 
dt — i — i — i 
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-w 
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X 
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» « 


where D. = 3x3 direction cosine matrix from body-fixed 
coordinates to inertial coordinates 

Equations 4.6) and 4.7) completely specify the attitude dynamics of Designs I 
and II. The inertia tensors of each design as well as a precise explanation 
of all coordinate frames and coordinate transformations are detailed in 


separate appendices. 
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4 . 3 Linearization of Equations of Motion 

When the direction cosines are perturbed slightly from their nominal 
values, it can be shown (Goldstein( 10) ) 


4.9) D . = D (I + E) E = 

— * — nom — — > — 


cos ijj - sin ijj 0 

sin cos 0 

0 0 1 


D 

— nom 


where D 

— nom 


-e 


z 

0 


-e 


-e 


x 

0 


4>= ) + ij/ , i)j 

o o 


If? rid/day 


1 = 3x3 


= nominal direction cosine matrix from body- fixed 
coordinates to inertial coordinates 
identity matrix 


(e ,e ,e ) = rotational errors about actual x,y and z body-fixed 
x y z 

axes, respectively with respect to each nominal 
(body-fixed) axis 

Using this, it is easy to show the actual angular velocity of the body- fixed frame 
with respect to an inertial frame, and its corresponding angular velocity, can be 
written as 


4.10) 

w 

w 

+ D 

» 

e 


— p 

~ — nom 

— nom 



4.11) 

« 

w 

» 

~ w 

+ w 

x D 

® T*1 

e + D e 

~p 

= — nom 

— nom 

— nom 

— — nom — 


where 


* 

- - 


■* ■ 


- -r 


e 


0 


0 


X 





£ = 

e 

w = 

0 

» 

w = 

0 


y 

— nom 

• 

— nom 



e 
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z 





■ — 


— 
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Sail torque can be expanded in a series, including only the nominal sail 
torque and a first order perturbation: 


4.12) N = N + D N. , 

— sails — nom — nom — body 

where N = nominal sail torque (hopefully zero, as will be shown) 

— nom ' 

N, = sail torque beyond nominal torque, in body-fixed axes 

— body 

Substituting 4.9) - 4.12) into the total equations of motion, including only 

terms to first order in linearized variables, and equating nominal trajectory 

variables to each other and perturbational variables to each other, the 

linearized equations of motion are found to be: 

I ( w xD e + D e) + 2(W I -I W ) D e 
—sat —nom — nom — — nom — —nom —pm —pm —nom — nom — 

4. 13) + (W* I . - I . W* ) D e+w xl D e+D ex (2 1 w 

— b — ab — ab — b —nom— —nom —sat — nom— — nom— —pm —nom 


1 , (W, + w- )) = D N, , 

— ab — b —nom — nom — body 


where I = I . + 21 

—sat — ab — p 



r 

ft 

- 







r « 


0 

-4, 

0 



0 

w 

az 

-w 1 

ay 


W 

ax 

w 

-nom 

•I*' 

0 

0 

wf = 
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and I and I are inertia tensors defined in an appendix, 
—pm — ab 

Two simplifying assumptions can now be made: 


i) I D *e » I (w x D e) + 2(W I - I W ) D e 
—sat —nom — —sat —nom — nom— —nom —pm —pm —nom — nom — 


-21 w xD e + w xl D e 
— p — nom — nom — —nom —sat — nom — 
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or, in words, the gyroscopic coupling of the satellite angular momentum 

vector with the angular momentum due to the annual rotation of the solar 

panels about their shafts, is neglected. Note each term on the right hand 

side of the inequality has a 365 day periodic component, while hopefully the 

term on the left hand side has most of its energy at much higher frequencies. 

ii) I D e » (W~* I . - I . Wf ) D e + D e x I ^ (wf + w ) 
—sat — nom — — b — ab — ab — b — nom — — nom — — ab — b — nom 

or in words, the gyroscopic coupling of the satellite angular momentum 
vector with the angular momentum due to the daily rotation of the antenna 
about its shaft, is neglected. Each term on the right hand side of the 
inequality has a 24 hour periodic component, while hopefully the term on the 
left hand side has most of its energy at much higher frequencies. 

Both statements can be summarized: it is assumed the sail torques 
can correct gyroscopic disturbance torques caused by antenna rotation about 
its shaft once a day.f-solar panel rotation about their shafts. These assump- 
tions will be checked later. 

The approximate equations of motion for the linearized variables are: 

4.14) I D e S D N, . 

—sat — nom — — nom — body 

4.4. Sensor Measurements 

4.4.1. Linearized Sensor Measurements 

Provided all pointing errors are small, earth sensors and sun 

sensors detect the center of the earth and sun, respectively, in body-fixed 

coordinates. Since £ is defined as the pointing error of actual body-fixed 

axes with respect to nominally sun-pointing coordinates, the components of 

e have the following physical significance 
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i) e -- rotation about the body-fixed x-axis , rota- 
x 

tion about the sun line 

ii) e -- rotation about the body-fixed y-axis, moving 
up-down with respect to the center of the sun 

iii) e -- rotation about the body-fixed z-axis, moving 
right-left with respect to the center of the 
sun 

assuming I e « 1 , e 1 , e « 1 . Sun sensors measure s e and e 

' x 1 y z 1 y z 

directly, but provide no information about e . 

Earth sensors must be used to find e indirectly. Consider the direc 

x 

tion cosines from the nominally earth pointing coordinate frame to antenna 
principal axes frame, D 


4.15) D 3 I + E' E’ = 
— a — — ) — 


z 

0 


-e' 


-e' 


where (e 1 ,e' ,e ! ) = rotation errors of actual antenna principal 
x y z 

axes with respect to nominal earth-pointing 


axes 


|e' J « 1 , | e ’ | « 1, j e ’ j « 1 is assumed. Using standard Euler 
x y z 


and 
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angles to go from antenna principal axes to body-fixed axes (see Appendix), 
it can be shown: 

e = e' (sin ib cost? + cos ibcos© sin*?) - e' sin9 sin<P 
x z y 

4.16) e = e' (sin lysin'? + cos tj* cos0 cos(| ) -e 1 sin© cost? 
y z y 

(V 4%/ 

e = e 1 cos ib sin0 - e' cosy 

Z Z y 

fJ A/ 

Since e' and e' are observable using earth sensors, they can be used to 

y 2 

calculate e ; as a bonus, earth sensor measurements provide a check on 

e and e measurements, and vice versa. These statements are valid 
y z 

provided the trigonometric weighting terms in 4.16) are not zero. 

Two effects complicate sensor measurements. First, both sensor 
noise and additive Gaussian amplifier noise corrupts the measurements. 
Second, since solar pressure dynamics are extremely slow, it is probably 
unnecessary to continuously monitor satellite attitude, but rather to sample 
it at a rate much faster than that of any disturbances . Equation 4.14) is a 
continuous time differential equation; in order to approximate it by a discrete 
time difference equation, the following state variables are defined 
4. 17a) e^kT) = e(kT) 

4.17b) e (kT) = e(kT) - e [ (k-l)T] 

4.17c) e (kT) = D _1 (kT) I (kT) D (kT) N , (kT) 

—3 — nom —sat — nom — body 

where k = 1,2,3,... 

T = sampling period 

and it is implicitly assumed the samples are taken at equally spaced time 


intervals. If it is assumed the sails exert no torque along the nominal 
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trajectory, then e^(kT) is approximately constant, and the difference equa- 
tion approximation to 4. 14) can be written as 



' £l (kTf 


"l_ I_ J_~ 


ej[(k-l) T] 

4.18) 

£ 2 (kT) 

= 

oil 


— 2 [ (k “ 1 ) T ] 


e 3 (kT) 


0 0 1 


— 3 1 (k _ 1) T] 


where C) = 3x3 all zero matrix 
I_ = 3x3 identity matrix 


The sensor noise is assumed zero (as in Chapter Three), while the amplifier 
noise adds to the observations: 


4.19) y(kT) = [I_ 0 0] 


e^kT) 

e 2 (kT) 

£ 3 (kT) 


+ n(kT) 


where n(kT) = samples of (amplifier noise) stationary Gaussian random 
process sample function 

}r(kT) = sensor measurements corrupted by amplifier noise 
If a minimum mean square error estimate of each of the state variables 
is desired, given noisy observations y(kT), a discrete time Kalman filter can 
be used. 

One approach to suboptimal estimation of £ (kT) and £ (kT) is to 

1 2 


assume the solar pressure dynamics, since they are extremely low frequency, 
will effectively low pass filter the noise; this is equivalent to ignoring the 
sensor noise entirely, and using for estimates of £ (kT) and £ (kT) equations 

A 4 

(4 . 17a , b) . Other examples of practical attitude estimation schemes can be 
found in MacLellan(5 9) Much(62), and Trudeau(68). 
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4.4.2. An Alternate Approach to Attitude Estimation 
One major difficulty in using sun and earth sensors is that neither the 

l 

roll -pitch -yaw angles, nor direction cosines or quaternions, are directly 
observable. However, if gyroscopes are used to measure angular velocity, 
then the direction cosines or quaternions are directly observable, and the 
results of Chapter HI can be applied. 

To see this, recall that if direction cosines are used to characterize 
orientation in space, then the differential equation which describes how the 
direction cosines evolve with time is 


d_ 

dt 


11 

d 12 

d 13 

21 

d 22 

d 23 

31 

d 32 

d 33 


-w 


w 


-w 


z 

0 


w 


-w 


w 


X 

0 


11 

d 12 

d 13 

21 

d 22 

d 23 

31 

d 32 

d,, 

33 

ntation in 

space , 


then the differential equation describing how the quaternion evolves with time is 

0 -w w 


d_ 

dt 


1 

2 


w 


-w 


-w 


z 

0 


-w 


w 


-w 


X 

0 


w 


w 


w 


X 


-w 


Since (w ,w ,w ) are directly observable using gyroscopes, these equations can 
x y z 

be integrated directly to provide all components of the direction cosine matrix or 
each component of the quaternion. Recall each component of the direction cosine 
matrix or quaternion was assumed observable in Chapter HI; since that is now the 
case, the attitude estimation method developed there can be used to estimate the 
roll -pitch -yaw angles, and these angles can be used directly to define errors 
with respect to a nominal trajectory, just as with’sun and earth sensors. 
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4.5. Solar Sail Torques 

An appendix derives the exact torque generated by each sail, under 
the assumptions of no sail shadowing, no radiation falling on the rear of the 
sails, and small attitude pointing errors. 

The force acting on a triangular sail such as in Designs I and II can 
be represented as a vector acting at a point two-thirds of the distance out 
the longitudinal sail axis, and this point is called the sail center of pressure 
(the distance from the sail hinge to the center of pressure will be called the 
center of pressure lever arm, or sail lever arm). 

Solar pressure has three distinct components: 

i) Reflected solar pressure, largest of the three, directed normal 
to the surface which sunlight strikes 

ii) Absorbed solar pressure, next largest of the three, directed 
along the sun line 

iii) Reradiated or thermal solar pressure, smallest of the three, 
directed normal to the surface which sunlight strikes 
Three simplifying assumptions are now made: 

i) Reflected solar pressure is the only contributor to sail torque 

ii) sin — b k cos b^ = 1 k = l,2,...6 

where b = sail k pitch angle =4pb ,b , or Ol (b = 15°) 

K lOO J O 

iii) (n « s) = (n • s) 3? (n. • s) k,j = 1,2, 6 

_ _ -k - -J - 

where n^ = sail k normal, in body-fixed coordinates 

^ = sun vector in body -fixed coordinates 

n = abbreviation for n, 

— — k 



92. 


or in other words, the projection of each sail's normal along the sun line is 
roughly the same. Using these assumptions, from the appendix it follows: 


N = A(1 -E ) (2 I/c) {n • s) 2 { (b -b + b -b, ) i cos (a) 

— body s — — 1346s 



» - 


m m 




0 


0 


L°J 

+ (b +b_+b >1 +b.) 1 sin (a) 
1 3 4 6 s 

1 

- (b +b ) 1 sin (a) 
ZDS 

0 

+ 



.0. 


1 

B» J 






sin 


V 

4.21) 

< b i + VVV r 

cos 

. 0 ■ 

+ (b +b,-b -b) 
3 6 14 

0 

Y _ 


where 


* — 

z + h - l cos a + 2r sin (a) 
s 

0 

- r cos a 

A = area of sail (all sails assumed identical 

E = emissivity of surface of sail 
s 

2 

I = sunlight pressure constant = 1 380 watts/meter 
c = speed of light in vacuum 

a = sail cant angle = 45° degrees 

j = distance from sail hinge on solar panel to center 
of pressure of sail 

r = radius of gyration of solar panel (see Appendix) 
h = distance from antenna center of mass to solar panel 
center of mass (see Appendix) 
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y = half length of solar panel (see Appendix) 
z = half width of solar panel (see Appendix) 
and all vectors are measured in body-fixed coordinates. Three cases of 
interest now arise: 



N* a 
— body 


o s 



1 


B — 
0 

cos(a) 

0 

-4b y 

0 



o 



0 

w * 


1 

u w 


+ 2 r sin (a) sin tp 


b) b = b = b = b, = b , b = b = 0 
1 3 4 6 o 2 5 


4.22b) N, b J = A(1 -E )(2I/c)(n . s) 2 {4b l sin(a) 
— body s — — os 


c) b = b = b = b = 0. b = b = -b 

1 3 4 6 2 5 o 


4.22c) N C , SA(1-E )(2I/c)(n • s) Z {2b f sin(a) 
— body s — — os 


•2b z 
o 


'o' 


* • 
1 

1 

+ 2 r sin(a) sin ij. 

0 

0 

B * 


0 

m> 4 



0 


1 

0 

+ 2 r sin(a) sin i|i 

0 

1 

J 


0 

» * 


If the sails are designed such that 

f cos(a ),Jt sin(a) » r sin a, y,z 
s s 

then x,y, and z torques can be generated in body-fixed axes, by cases 
a),b), and c) respectively. Sails 1,3,4 and 6 produce either x or y torques 
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at any instant, but not both simultaneously; sails 2 and 5 generate z torques. 

) 

If I is assumed diagonal, then 4.14) become 
—sat 


4.23) e = 


-1 2 - 1.2 
I cos ib + I sin ib 
xx yy 


2(1 * - I * ) sin ib cos 4* 
xx yy 


2(1 * -I ^ ) sin ib cos ib 0 

xx yy 


T -1 2 -1 2 
I sin ib+ I cos ib 
xx yy 


zz 


N 


where 


—sat 


xx 

0 I 


yy 

o 


0 

0 

[ 

zz 


(see Appendix for correct 
expression for I ) 

S3.U 


Note that if an x torque is generated in body-fixed axes, the off-diagonal 

entries in the gain matrix for N above will couple this torque into the 

- — body 

y axis. Cross coupling torques between x and y axes are unavoidable, 

because of the coordinate transformations between inertial and body-fixed 

axes. To minimize their effect, I and I must be chosen comparable to 

xx yy 

one another . 

4.6. Solar Sail Control Law 


At this point the analysis of the linearized equations of motion has 


body 


been reduced to a previously solved problem. Two standard control laws 
will now be discussed (for more details, see Athans (54) and Flugge-Lotz (56)). 
Cross coupling torques will be assumed negligible from here on. 



If it were possible to generate simultaneously three independent 
torques, the attitude error and its first derivative could be used in a control 
law 

4.24) N, , = N, a . sign(e +k e ) , sign(e +k e ) sign (e +k 

— body — body x x x — body y y y — body z z 


where 


sign (x) 


= ^+1 x — 0 

-1 x < 0 


Since x and y torques cannot simultaneously be generated, one possible 
modification to the above control law is 


4.25) 


If [e +ke 1^ je +ke | generate x torque 

If e +ke < e+ke generate y torque 

x xx yyy 


Finally, if a deadband is included in the control law for each axis, every 
"sign" function is replaced with a dbz 

f sign(x) | x I > db 


dbz(x) = 


db = deadband 


£ 0 | x | < db 

Equations 4.24) and 4.25) with deadbands completely specify a possible 
control law for either Designs I or II. The gains k ,k ,k allow the designer 
to trade ringing and overshoot in the system transient response for the 
period of a limit cycle when the satellite is in its nominal trajectory (see 
FlUgge-Lotz (56)). Another possible control law is the time optimal control 


law; a simple single axis control law is discussed in Athans (54). 
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4.7. Disturbance Torques 

Five sources of disturbance torques, in addition to the cross-axis 
coupling torques already mentioned, can drive Designs I and II away from 
nominal attitude. The largest disturbance is due to the earth's gravitational 
field trying to align the satellite's longest axis along a line pointing to the 
earth's center of mass. The next largest disturbance is caused by sunlight 
reflecting, absorbing, and reradiating from the antenna. Third largest is 
caused by the satellite center of mass being misaligned. Fourth is micro- 
meteoroid impact torque s . Finally, the smallest disturbance torques are 
due to the satellite interacting with the earth's magnetic field; these magnetic 
torques are assumed negligible. 

Gravitational torques on an arbitrary rigid body, or gravity gradient 

torques as they are referred to in the literature, are straightforward to 

calculate (see Nidney (65)). It can be shown 

0 
-I 

xz 
I 

xy 

where w^ = orbital frequency of satellite = (2 tt/ZA) radians /hour 

I ,1 = off diagonal elements in satellite inertia tensor 

xy xz 

computed in earth-pointing coordinates 
N® 

— gravity gradient = satellite gravity gradient torques in 


4.26) N 6 “ 3w 2 

— gravity gradient o 


earth-pointing coordinates 


I = satellite inertia tensor in earth -pointing coordinates 
—sat 


I e = A b I A* 

—sat — e —sat — b 


A* . A b T . A”' 1 
— b — e — e 


A = direction cosines from body -fixed axes to earth 


— e 


pointing axes in body-fixed axes . 
These torques can be written 


4 27) N — A 6 N C 

— gravity gradient — b — gravity gradient 

and are plotted for Design 1, in body-fixed axes, in Figure 4,3 under the 

assumption the satellite is moving in its nominal trajectory. 

Sunlight pressure torques due to the antenna-electronics body are 
calculated explicitly in an appendix. The reflected and absorbed sunlight 
pressure torques are plotted in body-fixed axes in Figures 4.4 and 4.5; 
more information must be known before the reradiated torque can be calculated 
explicitly. 

The calculation of center of mass misalignment torques, torques due 
to the center of mass of the satellite not coinciding with the center of mass 
of the antenna, is straightforward (Goddard (5 7)): 


4.28) N = R 


— cm — a /cm 


xF + M [ 2R x (w xR ) + R , x(w xR 

—sails sat — a/cm — p — cm — a/cm —a — « 


cm 


where 


N = center of mass misalignment torques in body-fixed axes 

— cm 

R . = vector from antenna center of mass to actual satellite 

— a/ cm 


center of mass 


Figure 4. 3 



( Figure 4.4 

Reflected radiation 



z - Torque 


Figures 4. 3-4. 5 

Nominal Disturbance Torques 
for Design I vs Time 


-0.0008 
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w = satellite angular acceleration in body-fixed axes 
“P 

w = antenna angular acceleration in body-fixed axes 
A typical magnitude for ^ a ^ crn one quarter inch (Goddard (56)). 
Since sail torque is proportional to the sail center of pressure lever arm. 


4 ‘ 29) ^sails a 1 S A,1 - E s ,,21/C > 


sin (a) 


cos (a) j 


The sails must be designed such that the sail torque is much greater than 
the center of mass misalignment torque, that is. 


f A(l-E ) (2 1/c) 
s s 


cos (a) 


sin (a) 


»M |R | [ 2 1 w ] + | w | ] 

sat — a/cm — p 1 —a 


f cos (a) » I R , I and 1 sin (a) » I R . I 
s '-a/cm 1 s '—a/cm' 

An analysis has been carried out to verify that if the sail cant and pitch 

angles are offset slightly from their nominal values, the resulting disturbance 

torques are negligible compared to the attitude correcting torques. 

Finally, the largest expected micrometeoroid torques appears to be 

two orders of magnitude less than the attitude correcting torques and hence 

is negligible (Goddard (57), NASA - Naumann(64' ) . 

4. 8 . Computer Simulation Results 

As a check on the foregoing analysis, a computer program was 

written to simulate the dynamics and kinematics of either Designs I or II. 

Figures 4.6 through 4.12 show the results of a simulation of Design I 


recovering from an initial pointing error. The initial error was due to the 
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satellite being rotated away from its nominal attitude by three successive 
rotations of ten degrees each, first about the z axis, thi n the y axis, and 
finally the x axis. Sun and earth sensor measurements are corrupted by 
additive Gaussian noise, and a sensor sampling rate of once each minute is 
used to determine the pointing error and its derivative (see 4.20)). The 
attitude control law described in 4.22) with a deadband is used to pitch the 
sails. Gravity gradient torques, antenna reflected and absorbed sunlight 
torques, and center of mass misalignment torques are included. Table 4.1 
summarizes the numbers used for these simulations. 

The simulation results lend credibility to the previous assumptions 
that gyroscopic coupling torques and cross-axis coupling torques are negligible 
compared to the sail attitude control torques. Note that for the example just 
discussed, the harmonic content of the error and its derivative is in the 
neighborhood of ten degrees per hour; this rate is much higher than the 
angular velocity of the solar panels about their shafts, but it is comparable 
to the angular velocity of the antenna about its shaft, roughly fifteen degrees 
per hour. Hence, there is some coupling between the angular momentum of 
the satellite and the angular momentum of the antenna, but the attitude control 
law is able to compensate for it. 

The results presented here are representative of over fifty different 
simulations of the attitude control equations for both Designs I and II . 

Satellite transient response to pointing errors such as those just described 
were checked at eight equally spaced times of year (starting with December 21), 
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Figures 4. 7-4. 9 

Computer Simulation Results of 
Design I- -Phase Plane Plots of 
Rotational Error Derivatives vs 
Rotational Errors 

Initial Errors are Due to Three 
Successive Counterclockwise 10° 
Rotations about the z-axis, then 
the y-axis, and finally the x-axis 
away from Nominal Orientation 
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Figure 4. 10 


Figure 4. 11 




Figures 4. 10-4. 12 

Computer Simulation Results 
of Design I--Solar Sail Torques 
vs Time 

Initial Errors are due to Three 
Successive Counterclockwise 10 
Rotations about the z-axis, then 
the y-axis, and finally the x-axis 
away from Nominal Orientation 


- 0.40 
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and at different initial times of day (midnight, six a.m. , noon, and six p.m.). 
Finally, most runs were made with the sails misaligned in cant and pitch, 
the satellite center of mass misaligned, gravity gradient disturbance torques 
present, antenna sunlight pressure disturbance torques present, and sensor 
measurements corrupted by Gaussian noise. These simulation results 
verified the results of a theoretical analysis: When the satellite design 
parameters are perturbed slightly from their nominal values, and when dis- 
turbance torques act on the satellite the attitude control torques are sufficient 
to keep the satellite aligned along its nominal orientation to within attitude 


sensor noise limitations. 
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Figure 4. 1.3 

Design 1 - - Simulation Parameters 

I. Inertia Tensor Components 

A. Antenna Principal Axis Inertia Tensor Components 

2 

I 300 kg-meter 

axx 

2 

I 225 kg-meter 

ayy 

I 225 kg-meter^ 

azz 

B. Solar Cell Body-fixed Axes Inertia Tensor Components 

2 

I 964 kg-meter 

xx 

2 

I 482 kg meter 

yy 

2 

I 482 kg-meter 

zz 

II. Sail Parameters 

2 

Sail area 3 meter 

1 - Center-of -pressure to sail hinge 5 meters 

E 0.05 (dimensionless) 
s 

III. Control Law Parameters 

x-axis deadband (same for all three axes) 0.1 degree 

K = K =30 (radians/min) * ; K =40 (radians/min) * 
x y z 

sensor noise (same for all three axes) 

mean 0 degrees standard deviation 0.03 degrees 
Cant angles 

a^ = 45degrees k = 1 , 2 , . . . 6 
Pitch angles 

= 15 degrees k = 1,2, ...6 
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IV. Dimensions 

Center of mass misalignment 

r , = (x . , y . , z . ) x / = y / 

a/ cm a/cm a/cm a/cm a/cm a/cm 

10 centimeters 

Antenna (distance from center of mass to surfaces) 


z 


a /cm 


r = 0.725 meters r = 0.5 meters x = 0.725 meters 
ax ay az 


Antenna surface areas 

2 

A, = A, = 2.05 meter 
1 o 

Solar panel 

r = 0.25 meter 
h = 1 . 0 meter 

V. Weight Breakdown (Goddard (57 
A. Antenna-electronics 


Antenna 

25 lbs 

Power control 

23 lbs 

Batteries 

50 lbs 

Electronics 

80 lbs 

Attitude control 

25 lbs 

Structure 

25 lbs 

Contingency 

19 lbs 


2 

= A^ = A^ = Aj. = 1.45 meters 

y = 0 . 25 meter 
z = 0.10 meter 
), M.I.T. (60)) 

B . Solar cells panels 


Bearing dxiive 15 lbs 

Solar cells 45 lbs 

Electronics 10 lbs 

Cables 8 lbs 

Structure 25 lbs 

Contingency 19 lbs 


Total 


247 lbs 


Total 


122 lbs 
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APPENDIX A 


Coordinate Frames and Coordinate Transformations 

Five different coordinate frames are useful in discussing the attitude 

control dynamics of Designs I and II. Each of the following four coordinate 

frames have their origin at the satellite center of mass: 

X -Y -Z --nominally earth-pointing coordinate frames, with 
e e e 

the X -axis aimed at the center of mass of the 
e 

earth, the Y^-axis perpendicular to X and in the 
orbital plane, and the Z^-axis parallel to the 
earth's axis of rotation 

X -Y -Z --the antenna -electronics principal axes frame, 

& SL sl 

nominally coincident with X -Y -Z 

e e e 

X -Y -Z --nominally sun-pointing coordinate frame, with the 
9 S S 

X -axis aimed at the center of the sun, the Y -axis 
s s 

perpendicular to X g and parallel to the ecliptic, and 
Z^-axis normal to the ecliptic 

X fe -Y b -Z b --the body-fixed axes frame, parallel to the solar cell 
panel principal axes frame and nominally coincident 
with X -Y -Z 


s s s 

The fifth coordinate frame has its origin at the center of mass of the kth 
solar panel (see Figures 1 and 2): 


X ^-Y^-Z -the principal axes frame for solar cell panel k, 

nominally parallel to X -Y -Z 

s s s 
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Direction cosines from one coordinate frame to another are differentiated by 

a combination of mnemonic subscripts and superscript?. For example, the 

direction cosines from X -Y -Z to X -Y. -Z. are denoted by 

a a a b b b 

c& 

D with the superscript "a" denoting antenna axes, the 

-Ij 

subscript "b" denoting body-fixed axes 
The direction cosines from earth -pointing coordinates to sun-pointing 
coordinates can be expressed in terms of Euler angles: 



COS 4j 

sin 4* 

o" 


"l 

0 

0 


cos ^ 

sin 

0 

D C = 
— s 

-sin i|j 

cos 4> 

0 


0 

cos 0 

- sin 0 


- sin ^ 

cosl^ 

0 


0 

0 

1 


0 

sin 0 

cos 0 


0 

0 

1 


where 


or 


(J£ = daily orbital position of satellite 
0 = angle of earth's axis of rotation with ecliptic 
ijj = annual orbital position of satellite 

ib = «[/ (t-t ) + ib ib= 2 t/( 365»24) radians/hour 
o o J 

0 = 23.5° 

(t _t o ) + = 2 jt/ 24 radians/hour 

t = initial time of day (midnight) 

= initial position in orbit (side of earth away from sun) 
i|i^ = initial time of year (December 21) 


where all initial times are nominal times taken for convenience in the com- 


puter simulation of Designs I and II, and easily changed to other values. 
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APPENDIX B 


Inertia Tensors 


In X -Y -Za the antenna-electronics inertia tensor is: 

2L 3 l 


axx 


i 

— aa 


ayy 

0 


0 


azz 


In body-fixed axes X, -Y, -Z, this becomes: 

b b b 


I = A* I A b 
— ab — b — aa — a 


.a A b -1 b T 
At = ^ = A 


— a 


In solar panel principal axes, X^-Y kth solar panel inertia 

tensor is: 


^Pk 


pxx 

0 

0 


0 

I 

pyy 

o 


o 

o 

i 

pzz 


.Oi body-fixed axes X^-Y^-Z^, for Design I, either solar panel inertia takes 


the form: 


I =1+1 

— p -pm — pk 


I = m D 1 I D 1 ” 1 
—pm p — nom —I — nom 


k - 


0 

-rh 0 


2 ,2 

r + h 


•rh 

0 

2 


m - solar cell panel mass 
P 
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h — i 

Figure B . 1 

Solar Panel Dimensions 
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D 


I 


— nom 


cos^ - sin ^ 0 

sin 'V sin^ 0 

0 0 1 

— 


In body-fixed axes for Design II either solar panel inertia takes the form: 


I = I + I , 

— p —pm — pk 


_ II _ II -1 

I = m D L D 

—pm p —nom -II — nom 


In 


h“ 
o h 


0 0 
2 0 
0 


n 

D 

— nom 


0 0 

cos"^ 0 ~sin^~J 
0 1 0 


-sm 


0 cos/ 


sin |( = sin^ - sin ^ 


Finally, the satellite inertia tensor becomes, in X, -Y, -Z, : 

b b b 


I =1+21 
—sat — ab — p 

This is graphed as a function of time for Design I, starting at midnight, 
December 21, and continuing for twenty four hours, on the following page. 
Numbers to do the graph are taken from Table 4.1. 


a (kg -met 


NOMINAL. DESIGN I INERTIA TENSOR 
IN SUN- POINTING COORDINATES 


— 

Ixx 

— 

i 1 

lyy 



1247 


947 


Time (hours) 

FIGURE B. 4 
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APPENDIX C 


Sail Torques 

Sunlight pressure striking a surface generates a force with three 

physically distinct components: 

F . = F + F + F 

— sunlight — reflected — absorbed — reradiated 

Only the reflected component will be dealt with here 


F 

—reflected 


A( 1-E ) 
s 


(21/ c) (n. 


s) 2 (-n) 


where 


A = surface area 
E = surface emissivity 

2 

I = sunlight pressure constants 1380 watts /meter 

g 

c = speed of light in vacuum = 3 x 10 meters/ sec 

^ = unit vector pointed at sun 

n = unit vector normal outward from surface 

The unit normal vectors for each sail in X, -Y, -Z, are: 

b b b 



cosa, 

k 

cosb, 

k 



cosa, 

k 

cosb, 

k 

= 

sina, 

k 

cosb, 

k 

k = 1,4 

^ = , 

- sina, 
k 

cosb, 

k 


L sinb k 




sinb, 

k M 








cosa 

cosb 2 



cosa 

5 

cosb,. 

—2 = 

sinb 2 


—5 = 


sinb,. 


sina 2 

sinb.. 



- sina^ 

cosb 

b 


^ ml 



— 




where 


a^ = cant angle with respect to x^ axis of sail k 
b^ = pitch angle about longitudinal axis of sail k 
The vectors from the satellite center of mass to the center of pressure of 
each sail are: 
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The torque generated by sail k is (in body- fixed axes): 

4 = -k*-k = A(1 - E 8> ,2I/c) ( £k‘ ^ ( 2k x £k> <» 

-body = -1 + \ + -3 + ^4 + -5 + -6 = total Sail tort * ue 


and where the dot product of the sun pointing unit vector with each sail normal 
is carried out assuming all vectors are expressed in body-fixed axes. 
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APPENDIX D 


Antenna Sunlight Disturbance Torqu es 
A. Reflected Sunlight Pressure Torques 

In antenna principal axes, vectors from the satellite center of mass 
to the surface of the antenna are: 
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The reflected sunlight force on each face of the antenna is 

sL = V‘- E ,k> (2I/ •’ ( - • ( -Sk> k = 1 6 

where 


\ ■ 
E ,k = 
2k f 


surface area of face k, k = 1 , . . . , 6 
emissivity of surface k 

unit normal vector outward from surface k in X -Y -Z 

a a a 


£ = sun-pointing vector in antenna principal axes 

The reflected sunlight torque due to each face is (in X -Y -Z ) 

2t cL 


4a = 2k x 4a ’ \ (1 - E sk> (2I/C) ( 2" Sk> % * 

The total reflected sunlight torque due to the antenna in body-fixed axes 

X -Y -Z^ is then 
b b b 


T 1 = A a T r total 
—ant — b —a 


T * , = T* + T* + T* +T T . + T* + T* 
— total —la —2a —3a —4a —5a —6a 
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B. Absorbed Sunlight Pressure Torques 

I 

Sunlight being absorbed by a surface it strikes causes a force to be 
generated. 

F , = A B (I/c) {s , n) (- s ) 

— absorbed — _ _ 

where 

A = surface area of surface 
B = absorptivity of surface 

2 

I = solar pressure constant = 1380 watts/meter 

g 

c = speed of light in vacuum =3x10 meters/sec 
j» = unit vector pointed at sun 
n = unit normal vector outward from surface 
The absorbed sunlight pressure force on each face of the antenna in antenna 
principal axes is then: 

-ka = A k B k (I/C) ( - # B) 

where 

B = absorptivity of surface k 

K 

The absorbed sunlight torque due to surface k in antenna axes is: 

W (I/C) 'i' V 

The total absorbed sunlight torque due to the antertna is the vector sum of the 


individual surface torques. 
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C. Reradiated Sunlight Pressure Torques 

Sunlight being reradiated by a surface it strikes generates a force: 


—re radiated ~ A (I R /c) ( ~- 


where 


A = Surface area of surface 


I_ = intensity of reradiated sunlight 
iv 


.8 


c = speed of light in vacuum = 3 x 10 meters/sec 
n = unit normal vector outward from surface 


To find the temperature on the surface must be found as a function of 


time, then the Stefan- Boltzmann formula must be used to find I 


R' 


Ir = E g t4 <‘» r 

= emissivity of 

^ = Stefan- Boltzmann constant 
T = surface temperature 


Final Note: When sunlight does not strike a surface, there is no 
reflected or absorbed force. Each of the forces in sections A and B for a 
given surface are multiplied by a function which is one when the surface is 
struck by sunlight and zero when it is not; this shadowing factor was included 
in the graphs of the reflected and absorbed antenna sunlight torques 
(Figures 4.4 and 4.5). 
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Chapter Five 

Summary and New Directions 

5.1. Summary 

The principal theme of this research program has been rigid body 
mechanics, as applied to space satellites. The study was broken down into 
three sections. 

The first section developed a four-dimensional space-time formulation 
of Newtonian mechanics. The equations of motion for a particle, a point with 
mass m, were found to be: 








-■ M « 0 

d* 

& ( My flofo) - 



where (x ,y , z ) are the spatial coordinates of the particle, (p ,p ,p ,m) is 
o oo xyz 

the linear momentum of the particle, and (F ,F ,F ) are the forces acting on 

xyz ° 

the particle, with all variables measured in a right-handed Cartesian inertial 
coordinate frame at the same instant of time t . The first four equations 
describe how the time rate of change of the linear momentum is related to 
forces, while the last six describe how the time rate of change of the total 
angular momentum of the particle, the moment about the space -time origin of 
the linear momentum, is related to forces. When no forces act on the particle, 
the four components of the linear momentum and the six components of the total 
angular momentum are all constant. Noether (21) showed through a Lagrangian 
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formulation of mechanics as well as a calculus of variations argument that 
these are the only ten quantities conserved when no forces act on the point 
mass; the formulation of mechanics pre sented here complements this view- 
point, and is perhaps more straightforward and offers greater physical insight 
into the nature of the conservation laws. 

The second section developed a new method for estimating the three 
angles which specify the spatial orientation of a rigid body, from observations 
of direction cosines and quaternions. The method presented is analogous to 
phase-locked-loop phase estimation, but is a generalization to three dimensions 
of phase-lock techniques, and is not simply three "one -dimensional" phase- 
locked loops. Two ways to implement this method, one based on direction 
cosines, the other on quaternions, were discussed for two specific examples, 
fixed orientation to space and constant rate of change of the three angles. The 
theoretical deterministic performance limits of each implementation were 
addressed: first, an analysis of the dynamics of each implementation for 
small estimation errors was carried out, and second, the steady state and 
stable equilibrium points of each scheme were discussed. When noise was 
included in each implementation the theoretical stochastic performance limits 
were quite difficult to pin down; in particular, Fokke r - Planck equations could 
be derived for each example under certain simplifying assumptions about the 
nature of the noise, but they could not be solved. However, for a special set 
of conditions, with a quaternion based method for estimating the unknown 
phases of the three angles as the angles change at a constant known rate with 
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time, the Fokker-Planck equation was solved explicitly for the steady state 
error probability density. All theoretical analyses, both deterministic and 
stochastic, were backed up by extensive computer simulation. 

The third section was concerned with the attitude control dynamics of 
two specific communication satellites, each using sunlight pressure to generate 
attitude control torques. Large reflecting surfaces, called solar sails, were 
attached to each satellite; by canting the sails in different directions, sunlight 
pressure torques were produced which corrected pointing errors. The equa- 
tions of motion for each design were derived, and then linearized about a 
nominal trajectory, and the linearized equations were analyzed using standard 
linear system theory techniques. The attitude of each design was determined 
from sun and earth sensor measurements. A heuristic bang-bang control law 
which governed when and how to cant the sails was developed. Five types of 
disturbance torques were identified, the largest of which is due to the earth's 
gravitational field, and were shown to be much smaller than sunlight pressure 
attitude control torques. Extensive computer simulation of the total equations 
of motion for each design, typically with a small initial pointing error, con- 
firmed the theoretical analysis of the linearized equations of motion. 

5.2. New Directions 

Several topics for further research are now outlined. These topics 
are grouped into the same three areas discussed in the main body of this 
report, in keeping with the central theme of developing a greater understanding 
into space satellite dynamics with applications to sunlight pressure attitude 


control. 
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1 . Deterministic Rigid Body Mechanics 

At the present time a mathematically rigorous theory of finite -dimen- ' 
sional linear systems exists, which has found wide application to many practi- 
cal problems (e.g. , problems arising in engineering and physics). At the 
same time, a mathematically rigorous theory of rigid body mechanics does 
not yet exist, but it would be useful in the design of practical attitude control 
systems for space satellites. Just as linear system theory has both qualita- 
tive and quantitative aspects, so would a theory for rigid body mechanics. 

A qualitative theory of rigid body mechanics would presumably draw on 
differential geometry and differentiable manifold theory. Lie algebras and Lie 
groups, as well as geometry and topology, in addressing such issues as con- 
trollability and observability of a space satellite, as well as questions of 
stability. One example in which such a theory might prove of immediate use 
is the design of control laws for despinning space satellites from high spin 
rates just after injection into earth orbit down to much lower spin rates. A 
second example where this type of theory might find immediate application is 
in the question of the stability of a dual- spin space satellite. 

A quantitative theory of rigid body mechanics would be concerned with 
at least two areas. The first is the application of digital computers and 
numerical techniques to questions such as numerical integration of rigid body 
equations of motion and optimal control laws for arbitrary space satellites. 

The second area is complex variable analysis. Elliptic functions and action 
angle variables play a key role in classical quantitative rigid body mechanics: 
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Hua(35) and Klein and Sommerfeld{l6) have indicated this role might be more 
significant than previously expected. Many problems in complex variable 
theory are analytically intractable, and perhaps can be solved numerically. 

Finally, how these threads of both qualitative and quantitative rigid 
body mechanics can be knotted together into a unified theory of rigid body 
mechanics, which would complement and lend insight to each other, remains 
to be seen. 

2. Stochastic Rigid Body Mechanics 

A natural extension of a deterministic theory of rigid body mechanics 
is to allow position, linear and angular momentum, and forces to be character- 
ized as random processes. One example of this was discussed in the section 
on attitude estimation, where it was observed that both direction cosines and 
quaternions could be considered as matrix exponentials of random processes. 
Three broad areas for further research into random processes and rigid body 
mechanics are now outlined: 

i) Several extensions of the work presented here, estimating 
the three angles which specify the orientation in space of a 
rigid body, come to mind. First, can the method presented 
here be extended to estimate the three angles if, for example 
the angular velocity is an arbitrary function of time, rather 
than the two specific examples discussed. Second, can that 
branch of information theory known as rate distortion theory 
be used to provide ultimate performance limitations for esti- 
mating the three angles in the two specific examples already 
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discussed, much as this theory has been used to point out 

I 

performance limitations in phase-locked-loop phase esti- 
mation (see Van Trees(47)). Third, how do quantization 
of observations (either direction cosines or quaternions) 
affect performance. Fourth, how is performance affected 
by sampling the observations at discrete instants of time, 
as compared to continuously observing spatial attitude, 
ii) A precise mathematical characterization of random 

processes found in rigid body mechanics is needed; since 
the problems encountered in mechanics are frequently 
nonlinear, mathematical rigor should aid in this charac- 
terization. Ito(36) and McKean(41) have carried out some 
preliminary work in this area. To be of practical use, 
their work must be related to observations made by actual 

attitude sensors; the section on solar pressure attitude con- 
« 

trol discussed sun and earth sensors, which are not yet 
adequately theoretically characterized, and which are not 
covered by Ito's or McKean's work. Another question of 
practical interest is how a force, which is a random 
process (e,g., disturbance torques in a gyroscope based 
attitude sensing scheme) influences attitude estimation 
(e.g., the three roll-pitch-yaw angles discussed earlier). 
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iii) A large body of mathematical literature exists on partial 
differential equations and properties of their solutions, 
both qualitative and quantitative . In particular , the 
Fokker-Planck equation has been extensively studied, as 
has abstract harmonic analysis, and the relationship of 
solutions of partial differential equations to representa- 
tions of those Lie groups and Lie algebras associated with 
spatial orientations. How this literature might help answer 
the question of what the optimal attitude estimation procedure 
is for processing noisy observations of space satellite orien- 
tation remains to be seen (see Evans(30)). 

3. Sunlight Pressure Attitude Control 

The intent of the third section of this report was to demonstrate the 
feasibility of a particular example of three-axis attitude control which used 
sunlight pressure to generate attitude control torques. Two specific designs 
were discussed, but perhaps there are other designs than those presented 
here which take better advantage of the technological constraints associated 
with space satellite design. The question of whether or not light weight rigid 
solar sails with the correct surface properties can actually be built has never 
been adequately answered; an adequate answer would be to actually construct 
and test in space a satellite which used torques generated by sunlight striking 
solar sails to control its attitude., and to the best of the author's knowledge, 
this has not been done for any earth orbiting space satellite. In particular. 
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one question that must be answered is how serious is sail flexing and bending 
apt to be and how much might this inference with attitud,- control dynamcis. 
Finally, can the same sails used for attitude control be used to station a 
synchronous equatorial orbit satellite over one spot on the earth's equator, 
in spite of irregularities in the earth's gravitational field. 


127. 

Bibliography 


A. Theoretical Rigid Body Mechanics 


1. R. Abraham, Foundations of Mechanics , W. A. Benjamin, N.Y. (1967) 

•9 — 

2. E. Bessel-Hagen, "Uber Die Erhaltungssatze Der Elektrodynamik," 

Math. Ann., 84, pp. 258-276 (1921) 

3. E. Cartan, "Sur Les Varietes A Connexion Affine et la Theorie de la 

Relativite Generalisee ," Ann. Ec. Norm. (3), 40, pp. 329-359 (1923) 

4. E. Cartan, Geometrie Des E spaces de Riemann , Gauthier- Villars, 

Paris (1946) 

5. H. B. G. Casimir, Rotation of a Rigid Body in Quantum Mechanics , 

J. B. Wolter's, the Hague, The Netherlands (1931) 

6. C. Che valley. Theory of Lie Groups I , Princeton Univ. Press, 

Princeton, N.J. (1946) 

7. H. Flanders, Differential Forms , Academic Press, N.Y. (1963) 

8. V. Fock, "Zur Theorie Des Wasserstoffsatoms," Zeit, Physik, 98 , 

pp. 145-154 (1936) 

9. H. Freudenthal, H. de Vries, Linear Lie Groups , Academic Press, 

N.Y. (1969) 

10. H. Goldstein, Classical Mechanics , Addison-Wesley , Reading, Mass. 

(1950) 

11. W. Greub, Multilinear Algebra , Springer- Verlag, N.Y. (1967) 

12. Y. Hagihara, Celestial Mechanics, Vol. I; Dynamical Principles and 

Transformation Theory , M.I.T. Press, Cambridge, Mass. (1970) 

13. S. Helgason, Differential Geometry and Symmetric Spaces , Academic 

Press, N.Y. (1962) 

14. N. Jacobson, Lie Algebras , Wiley Interscience, N.Y. (1962) 




128. 


15. F. Klein, Vorlesungen tjber Hohere Geometrie , Springer, Berhn (1962) 

16. F. Klein, A. Sommerfeld, Theorie Des Kreisels , Teuber, Leipzig (1897) 

17. S. MacLane , G. Birkhoff, Algebra , MacMillan, London (1967) 

18. G. McCarty, Topology , McGraw-Hill, N. Y. (1967) 

19. J. Moser, "Regularization of Kepler's Problem and the Averaging 

Method on a Manifold," Comm. Pure Appl. Math. 2^, pp. 609-636 
(1970) 

20. E. Nelson, Tensor Analysis , Math. Notes, Princeton Univ. Press, 

Princeton, N. J. (1967) 

21. E. Noether, "Invariante Variations Problems," Gott. Nachr., 

pp. 235-357 (1918) 

22. R. Penrose, "Structure of Space-Time," in Battelle Rencontre s , 1967 

Lectures in Mathematics and Physics , C. deWitt, J. Wheeler (Eds.), 
W. A. Benjamin, N. Y. pp. 121-235 (1968) 

23. A. Saenz, "On integrals of Motion of the Runge Type in Classical and 

Quantum Mechanics," U. of Michigan, Ph.D. Physics Dissertation 
(June, 1949) 

24. H. Samelson, "Topology of Lie Groups," Bull. Am. Math. Soc., 58 , 

pp. 2-37 (1952) 

25. H. Samelson, Notes on Lie Algebras , Van Nostrand Reinhold, N. Y. 

(1969) 

26. A. Trautman, "Noether Equations and Conservation Laws," Comm. 

Math. Phys., 6, pp. 248-261 (1967) 

27. F. Warner, Foundations of Differentiable Manifolds and Lie Groups 
Scott Foresman, Glenview, Illinois (1971( 

E. Whittaker, Analytical Dynamics of Particles and Rigid Bodies , 
Cambridge University Press, London (1965) 


28 . 



129. 


B. Attitude Estimation 


29. S. Bochner, "Positive Zonal Functions on Spheres," Proc. Nat. Acad. 

Sci. U.S.A., 40, pp. 1141-1147 (1954) 

30. J. Evans, "Chernoff Bounds on the Error Probability for the Detection 

of Non-Gaussian Signals," M.I.T. Elec, Eng. Sc.D. Dissertation 
(June, 1971) 

31. P. Frost, T. Kailath, "An Innovations Approach to Least-Squares 

Estimation- -Part III: Nonlinear Estimation in White Gaussian 
Noise," I.E.E.E . Trans. Auto. Con., AC-16, pp. 217-226 (1971) 

32. H. Furstenberg, "A Poisson Formula for Semi-Simple Lie Groups," 

Ann. Math., 77, pp. 335-386 (1963) 

33. R. Gangoli, "On the Construction of Certain Diffusions on a Differentiable 

Manifold," Z. Wahrscheinlichkeitstheorie , 2 , pp. 406-419 (1964) 

34. T. Hida, "Canonical Representations of Gaussian Processes and Their 

Applications," Mem. Col. Sci., Univ. Kyoto (A), 33 , 1> PP* 109-166 
(1960) 

35. I. K. Hua, Harmonic Analysis of Functions of Several Complex Variables 

in the Classical Domains , Trans. Math. Mon., Vol. 6, A.M.S., 
Providence, R. I. (1963) 

36. K. Ito, "Brownian Motion in a Lie Group," Proc. Japan Acad., 26 , 

pp. 4-10 (1950) 

37. K. Ito, H. McKean, Diffusion Processes and Their Sample Paths , 

pp. 269-276, Academic Press, N.Y. (1965) 

38. K. Ito, "The Brownian Motion and Tensor Fields on Riemannian 

Manifold," Proc. Int. Cong. Math., Djursholm, Sweden, pp. 536-539 

( 1962 ) 

39. A. Jazwinski, Stochastic Processes and Filtering Theory , Academic 

Press, N. Y. (1970) 

40. S. Karlin, J. McGregor, "Classical Diffusion Processes and Total 

Positivity ," J . Math. Anal. Appl., 1, pp. 163-183 (I960) 

41. H. McKean, "Brownian Motions on the 3-Dimensional Rotation Group," 

Mem. Col. Sci., Univ. Kyota (A), 33,1, pp. 25-38 (I960) 




130. 


42. H. McKean, "Brownian Motion With a Several -Dimensional Time," 

Theory Prob. Appl . , 8^,4, pp. 335-354 (1963) 

43. H. McKean, "The Bessel Motion and a Singular Integral Equations," 

Mem. Col. Sci., Univ. Kyoto (A), ^3,2, pp. 317-322 (I960) 

44. H. McKean, Stochastic Integrals , Academic Press, N.Y. (1969) 

45. F. Perrin, "Mouvement Brownien de Rotation," Ann. Ec. Norm. (3), 45 , 

pp. 1-51 (1928) 

46. D. Snyder, The State Variable Approach to Continuous Estimation , 

M.I.T. Press, Cambridge, Mass. (1969) 

47. H. Van Trees, Detection Estimation, and Modulation Theory, Part II: 

Nonlinear Modulation Theory , Wiley, N. Y. (1971) 

48. A. Viterbi, Principles of Coherent Communication, Me Gr aw Hill, N. Y. 

(1966) 

49. E. Wong, J. Thomas, "On Polynomial Expansions of Second Order 

Distributions," J.S.I.A.M., PP- 507-516 (1962) 

50. E. Wong, "The Construction of a Class of Stationary Markoff Processes," 

Proc, Symp. Appl. Math., Vol. 16, Stochastic Processes in 
Mathematical Physics and Engineering , A.M.S., Providence, R. I. 
(1964) 

51. E. Wong, Stochastic Processes in Information and Dynamical Systems , 

McGraw Hill , N. Y. (1971) ~~~ ~ — 

52. A. Yaglom, "Second Order Homogeneous Random Fields," Proc. Fourth 

Berkeley Symp. Math. Stat. Prob., Vol. 2, Contributions to 
Probability Theory , J. Neyman (ed.), U. Cal. Press, Berkeley (1961) 


r 



131. 


C. Solar Pressure Attitude Control 


53. J. Abel, "Attitude Control and Structural Response Interaction," JPL 

Report 32-1461 (1970) 

54. M. Athans, P. Falb, Optimal Control , McGraw Hill, N. Y. (1966) 

55. G. Fleischer, "Multi-Rigid Body Attitude Dynamics Simulation," JPL 

Report 32-1516 (1971) 

56. I. Flugge-Lotz, Discontinuous and Optimal Control , McGraw Hill, 

N. Y. (1968) 

57. Goddard Space Flight Center, "A Conceptual Study of a UHF-VHF 

Hybrid Satellite for a Delta Launch Vehicle," Internal Memo (1970) 

58. P. Likins, "Dynamics and Control of Flexible Space Vehicles," JPL 

Report 32-1329 (1970) 

59. D. C. MacLellan, H. A. MacDonald, P. Waldron, H. Sherman, 

"Lincoln Experimental Satellites 5 and 6," AIAA 3rd Com. Sat. 

Conf . , Los Angeles, California, 6-8 April 1970 

60. M.I.T. Center for Space Research, "Electrically-Steerable L-Band 

Satellite Antenna," Report 70-8 (1970) 

61. W. E. Morrow, ed., "Satellite Communications --Special Issues," 

Proc. IEEE, 59, 2 (1971) 

62. C. Much, N. Smith, F. Floyd, E. Swenson, "LES-8/9 Attitude Control 

System Numerical Simulation," M.I.T. Lincoln Lab. Tech. Note 
1971-10 (1971) 

63. National Aeronautics and Space Administration, Spacecraft Radiation 

Torques , NASA SP-8027 (1969) 

64. R. Naumann, "The Near Earth Meteoroid Environment," NASA TN D-3717 

(1966) 

65. R. Nidney, "Gravity Torque on a Satellite of Arbitrary Shape," ARS 

Journal, _30, 2, p. 203 (I960) 

66. L. Nguyen, "Three Axis Orientation and Stabilization of a Probe Using 

Solar Pressure," M.I.T. CSR Report T-70-4 (1970) 


132. 


67. R. Sohn, "Attitude Stabilization by Means of Solar Radiation Pressure," 

ARS Journal, 2_9, 5, p. 368 (1959) 

68. N. Trudeau, F. Sarles, B. Howland, "Visible Light Sensors for 

Circular Near Equatorial Orbits," AIAA 3rd Com. Sat. Conf . , Los 
Angeles, Cal., 6-8 April 1970 


133. 


Biographical Note 
Barton W. Stuck was born in 
He attended Grosse Pointe High School in Grosse Pointe. Michigan, graduating 
2nd out of a class of 781 in 1964. He received the degrees of Bachelor of 
Science and Master of Science, both in Electrical Engineering, from M.I.T. 
in 1969, and the Electrical Engineer degree from M.I.T. in 1970. During 
the 1968-69 school year he was supported in his studies by a Katherine Tuck 
Scholarship. From 1969 to 1972 he was a Teaching Assistant in the Electrical 
Engineering Department, except for the summer of 1971 and the spring term 
of 1972, when he held National Science Foundation Traineeships. 

His summer job experience includes work at the Bendix Research 
Laboratories in 1966, IBM Systems Development Division in 1967, M.I.T. 
Lincoln Laboratory in 1968, Hughes Aircraft Company in 1969, and M.I.T. 
Lincoln Laboratory in 1970. 



f 



